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SMALL BENDING AND STRETCHING OF SANDWICH-TYPE SHELLS 

By Eric Reissner 


SUMMARY 


A theory has been developed jor small bending and stretching 
of sandwich-type shells. This theory is an extension of the 
known theory of homogeneous thin elastic shells. It was found 
that two effects are important in the present problem, which are 
not normally of importance in the theory of curved shells: ( 1 ) 
The effect of transverse shear deformation and (2) the effect of 
transverse normal stress deformation. The first of these two 
effects has been known to be of importance in the theory of plates 
and beams. The second effect was found to occur in a manner 
which is typical for shells and has no counterpart in flat-plate 
theory. 

The general results of this report have been applied to the 
solution of problems concerning flat plates, circular rings, 
circular cylindrical shells, and spherical shells. In each case 
numerical examples have been given, illustrating the magnitude 
of the effects of transverse shear and normal stress deformation. 

The results of this investigation indicate the necessity of 
taking account of transverse shear and normal stress in sandwich- 
type shells, as soon as there is an order-of -magnitude difference 
between the elastic constants of the core layer and of the face 
layers of the composite shell. It was found that the changes due 
to transverse shear and normal stress deformation in the core 
may be so large as to be no mere corrections to the results of the 
theory without transverse core flexibility . 

The actual magnitude of the changes is greatly dependent on 
the geometry and loading condition of the structure under con- 
sideration so that no general rules may be given which indicate 
for which elastic modulus ratio the changes begin to be significant. 

Solutions of problems in the present theory may in general be 
obtained, by mathematical methods which are similar to those 
employed in the theory of plates and shells without the effect of 
transverse shear and normal stress deformation included. The 
present work does not include consideration of buckling and 
finite deflection effects. 

INTRODUCTION 

In this report an extension of the classical theory of small 
bending and stretching of thin elastic shells is considered. 
Instead of a homogeneous shell, investigation is made of a 
shell constructed hi three layers: A core layer of thickness h 
with elastic constants E c , G c , and v c and two face layers of 
thickness t with elastic constants E f , G f , and v f . In the 
developments certain restrictive assumptions are made which 
somewhat limit the general applicability of the results. In 
so doing formulas are obtained which are as compact as 


possible while still describing the essential characteristics of 
the sandwich- type shell. 

The thickness ratio tjh is assumed small compared with 
unity; at the same time the ratio E f t/E c h is assumed large 
compared with unity. This latter assumption means that 
the face material is so much stiffer than the core material 
that the contribution of the core layer to stress couples and 
tangential stress resultants of the composite shsll is negli- 
gible. It is known that for flat plates these assumptions 
necessitate the taking into account of the effect of transverse 
shear deformation. (See, for instance, reference 1.) The 
same would be expected to be true for curved shells, and the 
present report, therefore, gives a system of equations in 
which this effect is incorporated. 

A further effect which, it appears, has not been considered 
previously in the analysis of small deflections of sandwich 
structures is the effect of transverse normal stress deforma- 
tion. In the present report it is shown that this effect arises 
in a manner which is typical for shells and has no counterpart 
in plate theory. It may be likened, roughly, to what happens 
in the bending of curved tubes. 

The process by which the general results of this report are 
obtained is as follows : First, each of the face layers of thick- 
ness t is assumed to behave like a thin shell without bending 
stiffness. The loads applied to r these face shells, henceforth 
called face membranes, are of two kinds: (1) External loads 
and (2) loads caused by the stresses in the core layer. Next, 
the core layer of thickness h is assumed to behave like a three- 
dimensional elastic continuum in which those stresses which 
are parallel to the faces are negligible compared with the 
transverse shear and normal stresses. On the basis of these 
two assumptions three steps are carried out. First, the 
equilibrium equations of the core layer and of the face layers 
are obtained. Then an appropriate expression for the strain 
energy of the composite structure is derived. Finally, Castig- 
liano’s theorem of minimum complementary energy is used to 
obtain the relations which connect stress resultants and couples 
of the composite shell with the quantities which describe the 
state of deformation of the composite shell. 

The system of equations which is obtained in the foregoing 
manner is specialized for the following cases : 

(1) Flat plate 

(2) Circular ring 

(3) Circular cylindrical shell 

(4) Spherical shell with axisymmetrical deformation 
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In each case a number of problems are solved explicitly and 
the appreciable effect of transverse shear and/or normal stress 
deformation is illustrated numerically. 

This work was conducted at the Massachusetts Institute 
of Technology under the sponsorship and with the financial 
assistance of the National Advisory Committee for Aero- 
nautics. 

SYMBOLS 

h core-layer thickness 

t face-layer thickness 

£ 1; £ 2 curvilinear coordinates on middle surface of com- 

posite shell 

f distance coordinate measured along normal to 

middle surface of shell 

ai, a 2 coefficients of linear element on middle surface of 
shell 

Ri, R 2 principal radii of curvature of middle surface of 
shell 

N nmu direct stress resultants in upper face membrane; 

n— 1,2; m== 1,2 

N nm i direct stress resultants in lower face membrane 
fun, 'Pm tangential components of external load intensity 
on upper and lower membranes 
q u , <h normal components of external load intensity on 
upper and lower membranes 

ti f, r 2 f components of transverse shear stress in core layer 
o-f component of transverse normal stress in core layer 

Tniu, T nXi values of transverse shear stresses for f = A A/2 ; 

o-f „, ffn values of transverse normal stresses for f = ± A/2 

t nim values of transverse shear stresses at middle sur- 

face of shell 

Q u Q 2 transverse shear stress resultants 
N nm direct stress resultants parallel to middle surface 
for composite shell; 7i = 1,2; m=l,2 
Mnm stress couples for composite shell; 71 = 1,2; m=l,2 
p n tangential components of external load intensity 

for composite shell; n—l,2 

q normal component of external load intensity for 

composite shell 

s external load intensity term defined by equa- 

tion (22) 

II strain energy 

E f , G r , v elastic moduli of isotropic face-layer material; 

V= Pf 

E c , G c elastic moduli in transverse direction of core-layer 
material 

«i, u 2 effective tangential components of displacement 
of elements of composite shell 

w effective normal component of displacement of 

elements of composite shell 

di, /3 2 effective components of change of slope of normal 
to middle surface of composite shell 
component of strain (e fm = o- ! - m /£' c ) 

C*-—2tE f 

D* = (l/2)t(h+t) 2 E f 

C—C*/{1 — r 2 ) 

D bending stiffness factor (D = D*j(l — r 2 )) 

x, y Cartesian coordinates in plane of flat plate 
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r, d polar coordinates in plane of flat plate 

a radius of circular ring, cylindrical shell, and spher- 

ical shell 

x, 6 surface coordinates on cylindrical shell 

X], X 2 , X 12 parameters defined by equation (63) 

H = ir/l 

l half wave length of sinusoidal load distribution 

mi, m 2 quantities defined by equation (197) 
k complex quantity defined by equation (200) 

< p , 6 surface coordinates on spherical shell 

co quantity defined by equation (74) 

K parameter defined by equation (190) 

I— GENERAL THEORY 

STATICS OF SANDWICH-TYPE SHELL 

In order to derive a complete system of equations for the 
shell composed of face layers and core layers it is necessary 
first to consider separately the statics of the face layers and 
that of the core layer of the shell. Combination of the results 
obtained for these two components of the composite structure 
must and will lead to those differential equations of equilib- 
rium which hold for elements of a shell, whether this shell 
is of homogeneous or nonhomogeneous construction. In 
addition, however, relations are obtained which are charac- 
teristic of the sandwich-type shell. 

Coordinate system on shell. — A curvilinear coordinate 
system (£ 1; £ 2 , f) is introduced as follows: Let & and £ 2 be 
coordinates on the middle surface of the composite shell and 
let f be the distance of a point of the shell from its middle 
surface, measured along the normal to the middle surface. 
In order that this system of coordinates be an orthogonal 
system, choose the &, £ 2 curves as lines of curvature on the 
middle surface (in the case of shells of revolution the lines 
of curvature are identical with the meridians and parallels 
on the middle surface) . 

The linear element in the forgoing system of coordinates 
is of the form 

ds 2 =m 2 (l +0 dtf+a 2 2 (l W+df 2 (1) 

where a, and a 2 are the coefficients of the linear element on 
the middle surface and 7?, and R 2 are the principal radii of 
curvature of the middle surface (see fig. 1). Formulas for 
the calculation of the quantities a n and R n are contained in 
texts on differentia] geometry. They are collected, together 
with other results, in reference 2, which deals with the 
theory of homogeneous thin shells. 

Statics of face layers. — The face layers are treated as 
thin shells of thickness t and it is assumed that the bending 
stiffness of these thin shells about their own middle surface 
may be neglected. (This, of course, means that no local 
buckling phenomena are considered in the present work.) 
Because of this neglect from now on they will be designated 
as face membranes. 

The middle surfaces of the face membranes evidently are 
given with reference to the three-dimensional system of 

curvilinear coordinates by f=i (A-f-f) and f = — i (A-j-t). 

Li Li 
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Figure 1. — Element of composite shell, showing coordinates and dimensions. 


From equation (1) it follows that the linear element on the 
middle surfaces of the face membranes is given by 

*■=«.■(' ± $)’ <*«.■+.,■ 0 ± h M) iif (2) 

The components of external load intensity on the upper 
and lower membranes are designated by p lu , p 2u , and q u and 
by p u , Pih and q h respectively (fig. 2). The core-layer 
stresses which act on the upper and lower membranes are 
given as tk„, and <rj- M and by t 2 (i, and <tj- ^ respectively. 

Finally, the direct stress resultants in the upper and lower 
face membranes are designated by N n u, N 12u , N 21u , and 
N 22u and by N 1U , Nm, N iU , and N 22l , respectively (fig. 2). 



Figure 2.— Element of composite shell, showing location and orientation of stress resultants 
in face layers and core layer and orientation of external loads. 


There are then three equations of force equilibrium for the 
elements of each of the two membranes. Writing 



the equations for the upper-face membrane are the following: 1 


ba 2u N U u i do;iu-^2iu 


aiu<X2u{Plu — Tifa) =0 


AT dotlu AT b(X 2u 

IV 12 u IV 221 / 


5*1 


da 2 u N 12u t b<x lu N 22u 


-N 21 U 


OL\u0^2u 


af 1 1 &?i 

«1K «2„ (p 2u — r 2 f „) = 0 
Nnu I Nm u 

j-( i+ ys!) z-O+U!) 


bot 2u at u 

IVll u 




-1u~r 0{ U 


(4) 


(5) 


= 0 (6) 


The corresponding equations for the lower-face membrane 
are 


cioi 2 iNi\i 

2>*i 


t>auN nl 

5^2 


FA 12 ; 


daw 

£>£2 


N 22 i 


ba 2 i . 
2>€i 


«u« 2 i (pu+i in) — 0 

(7) 


ba 2 iNi 2 i , bauN 22 i , 


d£i 


d£ 2 


^y + « u «2! (p 2 f + T 2 n) = 0 

(8) 


ana« 


N u 


Nv 






= 0 (9) 


As bending moments and transverse shears are assumed 
not to be acting in the individual membranes the moment 
equilibrium equations become the symmetry relations 


Ni 2u —N 2 iu 

N l2 ,=N 2 u 

Before analyzing the state of stress in the core layer it is 
convenient to see what relations follow from equations (4) 
to (9) for the composite shell. 

Statics of composite shell. — It may be seen that, in view 
of the fact that all face-parallel stresses in the core layer are 
neglected, the following expressions for the face-parallel 
stress resultants and couples of the composite shell are 
obtained: 


]■ (10) 

j 




i These are obtained from the. corresponding equations of reference 2 with a n changed to 
a n% and with stress couples and transverse shear stress resultants omitted. To make up for 
this omission, the loads on the two membranes arc assumed to act at their middle surfaces, 
this means terms of the order t/R are neglected (but not terms of order h/R) . 
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N 2 i=( 


12a + ^1 


12 1 

(13) 


N 2 = ( 

'l+‘- 

m) N 

22a + (^l 


22 1 

(14) 

Mu 

h-\-t 

2 

[(* 

, A>+A 
1 2 rJ 

Nuu — 

O-wO 

M,«] (15) 

M 12 

h-\~t 

2 

[(■ 

A+f\ 
1 2 RJ 

Ni2u — 

0 2ft) 

A7i21 

(16) 

M n 

h -\~t 

2 

K 1 

+ h+t ) 
1 2 RJ 

N 12 a 

0-&0 

N 121 

(17) 

m 22 

_h + t 
2 

[0 

A+<\ 
1 2 RJ 

^N 22 a — 

0 \t) 

a 22( ] 

(18) 


In the same way the following expressions are obtained 
for components of external force and moment intensity: 



f 2E 2 )( 1 + 2I?;)^“ + ( 1 

h+t )( 1 . 

2 rJv 

2Rj Pnl 



(19) 

►wDO+wV-K 1 - 

-2r)(‘ 

h T t\ 



(20) 


(-$)(' 

A-(-f\ ~i 
2R l P nl \ 


( 21 ) 


Further, a load term of the following form will be 
encountered: 








which bears a relation to equation (20) similar to that which 
equation (21) boars to equation (19). This last term would 
represent, for a homogeneous shell, the average transverse 
normal stress at any station of the shell, assuming that the 
loads q u and q t alone are responsible for this stress. For a 
homogeneous isotropic shell this term is of no importance. 
For a sandwich-type shell, as will be seen, it may sometimes 
be of importance. 

In order to obtain force and moment equilibrium equations 
for the composite shell the face-membrane equilibrium 
equations (4) to (9) are combined suitably. Adding equa- 
tions (4) and (7), and (5) and (8), respectively, the two equi- 
librium equations for the force components parallel to the 
middle surface of the shell are obtained. In order to reduce 
them to known form (see reference 2) the following relations 
are used between the core-layer-surface shear stresses r n 
and T ntt and the transverse shear stress resultants Qi and Q 2 . 



(23) 






(24) 


Equations (23) and (24) will subsequently be shown to be 
in agreement with the usual definition for the transverse 
shear stress resultants by consideration of the stress distri- 
bution of the core layer. 

With equations (23) and (24), there are obtained by com- 
bination of equations (4) and (7), and (5) and (8) — carrying 
out addition as well as subtraction — the following four 
equations : 


Pi ^ 

r 

b N12 

£»^2 

ba 2 Ni2 , 

bcuNp, Ar 

ba 2 

Pi H 

P2 

r4.v2i 

Pi 

ba 2 M n 

. C)aiM 2 i 

1 . 71 /f 


Pi 

‘ P2 

^ 12 p2 


|_ ^OL\J\d^22 

f 



12 l7r N22 Wi + 011012 (r 1 +Pi )~ 

2 ( J +252 ) = 


-Nn 5^ + aia 2 


2 

22 Pi 

daj 


+ «i <x2(mi — Qi) — 0 


-M u £j^ + C*l“ 2 (W'l 2 — Q 2 ) 0 


(25)* 


(26)* 

(27) * 

(28) * 


Two further equations are obtained by adding and sub- 
tracting, respectively, equations (6) and (9). Adding equa- 
tions (6) and (9) and taking account of equations (11), (14), 
and (20), there follows: 


«1«2 


'/Nil ,Nn 
A Ri R* 



1 + 


h d~ t \ / . , A-f A 
2RJV^2RJ 


Otu — 



(29) 


In order that this reduces to the correct equation of trans- 
verse force equilibrium as given in reference 2, one must have 


/ A + A/ h+t\ "] 1 dcnfe 

V 1 “2^A 1 “2i7;r r 'J = ^r + ^r 


h -\-t\ 

" 2R V A f “ 


(30) 


Equation (30), just as equations (23) and (24), can again be 
verified independently by consideration of the state of stress 
in the core layer. On the basis of equation (30), equation 
(29) is written in the form 


d<*2<2i , dcnQs 


d\ do 


/AT 

2 \ Ri 


11 i N 22 


4~ g- ) + aia 2 g = 0 


(31)* 


The last equation, use of which is required for the sand- 
wich-type shell and which has not previously been given, 
is obtained by subtracting equation (9) from equation (6). 
Taking account of equations (15), (18), and (22), there 
results 


2a.\CL2 

K-\-t 



^)-2 a ia2 s + 





0 




(32) 
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Provisionally, there is written 


0+^)0 






)0 


h-\~t 
'2 R 


^ <Tfi — 2 

(33) 


and it will subsequently be shown that represents the 
value of erf at the'"middle'surface of the shell. Combining 
equations (33) and (32) yields 


a ^ + h+t 



(34)* 


Equation (34) has no relation to the sixth equation of 
equilibrium for an element of the shell which expresses the 
condition of moment equilibrium about the normal to the 
middle surface. That equation which, as is known, is an 
identity when resultants and couples are expressed in terms 
of stresses does not occur in the present derivations, or 
rather it is contained in equations (12), (13), (16), and (17), 
which give explicitly the slight differences between N l2 and 
Nn, and M. VI and di- 
stress distribution in core layer. — In order to verify 
independently equations (23), (24), and (30), as well as for 
the subsequent derivation of appropriate stress-strain rela- 
tions, it is necessary to determine the distribution of stress 
in the core layer. 

Assuming that the components of stress <r u <r 2 , and r 12 in 
the core which would contribute to stress resultants and 
couples of the composite shell are of negligible importance, 2 
these components of stress may be set equal to zero and 
only the components of transverse shear stress and trans- 
verse normal stress r lf , r 2 f, and <r f may be retained. The 
differential equations of equilibrium for these three remain- 
ing components of stress in the system of curvilinear coordi- 
nates defined by equation (1) are obtained, from the general 
form of these differential equations in reference 3, in the 
following form: 


O+AX^A)^- 

f f d / airi jm \ . d / a x T 2t m \~| 
“ia 2 + f/^i/ + 


(40) 


The transverse shear stress resultants Q x and Q? are 
obtained from equations (38) and (39) in the form 


X (.h+t)l 2 
-( h+t)/2 


Tni 




(h T i) 

T n$m 



(41) 


The integration must be extended over the thickness of the 
core layer and also over half the thickness of the face layers, 
in accordance with the prior assumption that the stresses 
T n [u, t „ f( , o-f,,, and o-n may be taken to act at the middle 
surfaces of the respective face membranes. 

Now, as intended, the proof is carried out of equations 
(23), (24), and (30), which were used to obtain the differ- 
ential equations for the composite shell. 

To verify equation (23), from equations (38) and (39) 
for the left-hand side of equation (23), the following equa- 
tion is obtained : 


r n$m T nX m ^ 1 ^ T n£m 

, , h-\-t h-\-t R n 1 / fc-j-f V 

+ 2 R n 2 R n \2 R n ) 

and this, in conjunction with equation (41), verifies equation 
(23). 

To verify equation (24) in the same manner, from equa- 
tions (38) and (39) for the left side of equation (24), the 
following equation is obtained : 


d^ 

df 


d 

dr 

[O+^O+AXAH 

(35) 

d 

ar 

K 1+ A) ( 1+ A) T2r ] =o 

(36) 


;)( 1+ £) ff 0 + M“ 2 ( ]+ A) 

1 7 lf J + 


&[ a >( l+ £) T2f ] =o 

(37) 


The values of the three stress components at the middle 
surface (f=0) are designated by the subscript m. Integra- 
tion of equations (35) to (37) then gives 


T\( = 


T 2 r= 


. 7 l lm 

(i + f/f^i) 2 (i + f/fy 

m 

(1 + f/E 2 ) 2 (1 + f/i?i) 


(38) 

(39) 


2 It is for this purpose that the order-of-magnitude relation hE c /tEj<& 1 is assumed. 


h-\-t , 

( Tntm | 

T n£m \ 

(A + i) t 

2 ( 

l+~ 
V + 2 R n 

1 h-\-t 

2 Rj 



and this, in conjunction with equation (41), verifies equa- 
tion (24). 

To verify equation (30), equation (40) is used to write for 
the left side of equation (30) 


d / a 2 r^ m \ d / a X T 2{ m \ . d / a 2 rif m 


v 



f OtlT# m \ 


2a 2 T l{m ~ 


2ai7 2 j- m 


h-^-t 

\ 2R 2 J 

| ^ 

L-csgyJ 

i_/A+«v 
L \2i?2/ _J 



and this, in conjunction with equation (41), verifies equa- 
tion (30). 

The section on the stress distribution in the core layer is 
concluded by listing the form which equations (38) to (40) 
for the stresses in the core layers assume for “thin” shells, 
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that is, for shells for which h/R< Cl. From equations (38) 
and (39), in conjunction with equation (41), it follows that 


_ 

U( -h+i 


y 



(42)* 


From equation (40), in conjunction with equation (41), it 
follows that 


a 1 a 2 <r t =a 1 a 2 <T tm — 


) 


(43a) 


It is necessary to note for some of the following considera- 
tions that, in view of equation (31), instead of equation (43a) 
there may be written 

f (Nw | N 22 

o-i— <rr™ h+t\Ri R 2 V (43b) 

It is seen that in this approximation the transverse shear 
stresses are uniform across the thickness of the core layer, 
while the transverse normal stress is composed of two terms, 
one uniform across the thickness and the other varying 
linearly across the thickness. 

No further calculations are needed with reference to the 
state of stress in the composite shell. The next step is to 
complete the system of differential equations for stress re- 
sultants and couples by deriving an appropriate system of 
stress-strain relations. 

STRAIN ENERGY OF SANDWICH-TYPE SHELL 

In calculating the strain energy of face membranes and 
core layer it is assumed that both are isotropic and elastic, 
with elastic constants E f , v f =v, G f —E r l 2(l + r) and E c , v c , 
G c = EJ2 (1 -j- v c ) . Poisson’s ratio for the face membranes is 
written without a subscript, because, in view of the assumed 
stress distribution, there is no explicit occurrence of Poisson’s 
ratio v c for the core layer. 

The strain energy for the composite shell is the sum of the 
strain energies for the face membranes and for the core layer 

n=n r +n c (44) 

For the purpose of obtaining stress-strain relations, both 
11/ and II C are expressed in terms of stresses rather than in 
terms of strains. 

Strain energy of face layers. — Considering that the element 
of area on the middle surfaces of the membranes is of the 


form aia 2 (l± 7 ^^) (l and that the stresses 

in the membranes are the stress resultants divided by the 
membrane thickness t, there is obtained, from well-known 
principles, the following relation: 

n /=i J JJ=r [NuJ+N 22 J- 2vNnuN UH + 2(1 + p)N 12 J] X 

if J tWf i N "‘ 2 + N ™- 2vNniN ™ +2 ( 1 + v)N '* 2] x 

( i -&0( i -&0 W6,i6 (45> 


Equation (45) is transformed into an expression containing 
stress resultants and couples of the composite shell by means 
of equations (11) to (18) which lead to the relations 



with corresponding formulas for N ]2 and N 22 . Note that 
equations (46) and corresponding equations can be used to 
calculate the stresses in the two different face membranes, 
once stress resultants and couples in the composite shell are 
known. 

In what follows attention will be restricted to cases in which 
h/R<g. 1. Then, with the two constants C* and D* defined by 


C* — 2tE t 
D*=^t(h+t) 2 E f 


(47) 


the following expression for 11/ is obtained : 

II/=iJJji [N n 2 +N 22 2 — 2 vN n N 22 + 2(1 + p)fVi 2 2 ] + 

i[M u 2 +M 22 2 - 2rM„M 22 +2(l+r)M 12 2 ]| (48)* 


It may be remarked that equation (48) could have been 
given directly, by analogy with known results for the 
isotropic homogeneous shell. 

Strain energy of core layer. — With the stresses a l; <j 2 , and 
T ,2 assumed to vanish, there results for the strain energy of 
the core layer 


k+t 




Q't 2 1 


G c 


) ( i+ £) ( i+ £) 


(49) 


Again the terms f/I? compared with unity are neglected and, 
consistent with this neglect, the values of the stresses r n{ and 
<r f are taken from equations (42) and (43). 

The value of <r f may be chosen from either equation (43a) 
or equation (43b). The form of the results depends some- 


what on which of the two equations is selected, in the sense 
that the meaning of the deformation quantities which are 
to be determined depends on which of the two equations is 
taken. This question is decided in the following manner: 
As all resultants and couples enter the expression for the 
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strain energy only as themselves and not in differentiated 
form, except when equation (43a) is used, the selection of equa- 
tion (43b) for <T f is proposed, thereby excluding derivatives of 


stress resultants and couples from the expression for the 
strain energy II. 

Introducing then equation (43b) into equation (49) yields 


J J 0Sr + % 2 ~ 2 )] w&rffii 


h -|- 1 \ R\ 1 R 2 

The. integration with respect to f is carried out and equation (50) becomes 


"•4 


r n< 2 i a +< 2 2 2 

h-\-t 

_ 2 _i , Af 22 vni 

j J Uh+t)6 c 

h E c 



(50) 


(51)- 


It was to be expected that the terms containing the 
modulus of rigidity 6 C would occur in the foregoing form. 
The contribution of the present report up to this point, 
besides giving the new equation (34) for cr Sm , is thought to 
he the determination of the form in which the effect of 
transverse normal stress deformability manifests itself in 
the strain energy of the sandwich shell. 

STRESS-STRAIN RELATIONS FOR COMPOSITE SHELL 

In what follows a system of stress-strain relations for the 
composite shell is obtained by the use of Castigliano’s theorem 
of minimum complementary energy. The manner in which 
the theorem is used here appears to have been employed 
first by Trefftz (reference 4) for the purpose of avoiding 
geometrical considerations in the derivation of the stress- 
strain relations for thin homogeneous shells with small 
deformations, without consideration of the effects of trans- 
verse shear and normal stress deformation. 

Assuming for the present purpose that all boundary 
conditions for the shell under consideration are stress con- 
ditions, the theorem consists in the statement that among 
all statically correct states of stress the actually occurring 
state of stress makes the strain energy of the system a 
minimum. In the application of the theorem the fact is 
taken into account that statically correct states of stress 
only are to be compared, by means of the Lagrangian 


multiplier method, before minimizing II an integral is 
added to it which contains the six equilibrium equations 
(25) to (28), (31), and (34), each of the six equations multi- 
plied by a Lagrangian multiplier. It can then be shown, by 
using Castigliano’s theorem with prescribed boundary 
displacements instead of with prescribed boundary stresses, 
that each of the six multipliers has the meaning of one of 
the displacement quantities which occur in the shell problem. 3 

With the foregoing understanding of the meaning of the 
multipliers, the multiplier of equation (25) is designated by 
u,; that of equation (26), by u 2 ; that of equation (27), by 
0i; that of equation (28), by 0 2 ; that of equation (31), by w; 
and finally that of equation (34), by k. It is known that 
u lt u 2 , and w represent the effective components of displace- 
ment in the £,, f 2 , and f directions, respectively. Further, 
it is known that 0, and 0 2 represent the angles through which 
the normal to the middle surface of the shell turns toward 
the £i and f 2 curves, respectively. There is no immediate 
simple geometrical interpretation for k and, while such inter- 
pretation in terms of an average transverse normal strain 
might be deduced herein, k is considered as an auxiliary 
variable presently to be eliminated. 

Combining now equations (44), (48), (51), and (25) to (28), 
(31), and (34) in the manner indicated, the following varia- 
tional equation results: 


5 If/ ^*+Nn t -2vN ll N»+2(l + p)N u t \ + 

~jy* [M u 2j t M . 22 — 2rM 1 iM 2 2+2(l+r)M 1 2 2 l + 


Qi 2 +<2 2 2 , 

, h-\-t 

(h-^rtyGc 

. 


OTm 2 + 


_ 

12 Vi?! ~ r i? : 


u i N22 
+ ~e> 2 


aia 2 d^\d^2~\~ 


[ 
[ 
[ 
[ 


ba 2 Nn , bctiNt 


a«i 


d£ 2 


21 , AT Cl^l AT C)d2 . f vs, 

biVi 2 577 77^-t-ai<* 2 ( -p- 


U. 


Pi 


. 5{i 

f H 

hi 


do!iM 2l 

1 

. 


“T 

5q!2-^12 

1 ba l M 22 

1 

. <>«. 

b£ 2 

t- 

ba 2 Qi . 

5«iQ 2 

-A . C 



d£ 2 

AT 
Nn 




+ 

+ 


(t +3 ’ 1 )] 

_2_iVii §S + “ i “ 2 (i +P2 )] 

M 12 1| 1-M 22 |^+« 1 «,(m.-<2 1 )]+ 

.M21 — Mn £jjr + ai« 2 (w 2 — Q 2 )J + 

/Nu,N» VI, 

ni?7 + i^- 2 ;J + 


k ( R [ l + £*") jd&dfe— 0 


(52) 


3 For the special case of the flat plate this has been carried out explicitly in reference 1. For the case of the homogeneous shell, without effect of transverse shear and normal stress defor- 
mation, the proof has been given in reference 4. The proof for the more general case which is here considered is not included as it does not offer any clearer insight into the problem and tends 
to lengthen the analytical discussion. 


888244—60 2 
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The variations in equation (52) are carried out and deriva- 
tives of variations in the double integral are eliminated by 
integration by parts. The line integrals along the boundary 
which occur due to this integration by parts vanish, because 


it has been assumed that all stresses are prescribed at the 
boundary and therefore their variations vanish at the bound- 
ary. 

The resultant variational equation is 


j SN 11 

<N 

£ 

& 

1 

£ 

1 

1 dui 

U 2 CiCKi 

w 

h+t j 

Wn , iV* 

L c* 

«i d£i 

ai« 2 d£ 2 

'W 

'12 E c Ri ' 

8 N 22 

r N 22 vN u 

1 du 2 

Ui d<* 2 

w 

h+t j 

Wu 


L c* 

do d£ 2 

ai<x 2 d£i 

R 2 

‘ 12 E c R 2 ' 

V if 1 


SNi 2 

f2(l -j- v)N \ 2 

1 bui 

1 U, 

d«i 

1 

d u 2 | u 2 

da 2 "l 

+ 

L c* 

a 2 d£ 2 

' di a 2 

^ ^2 

di 

d£i 1 0 : 10:2 

dfj 

5 Mn 

■A^ll vl \^22 

j_ M> 

02 

dai 

1 

1 k ‘ 

]+ 


D* 

d£i 

«1«2 

d£ 2 

(h-\-()Rl Ot\Ot 2 _ 


dM 22 

r~. a/ 22 — vi\£\ 1 


01 

da 2 j 

, 

1 k - 

]+ 


_ D* 

«2 d£ 2 

«1«2 

d$. 

(h -{- 1 ) R 2 


8 M 1 2 

2(1 -f- v)Mi 2 

1 Ml 

1 01 

bai 

1 

d/?2 | 02 

da 2 

+ 

D* 

a 2 d£ 2 

dld 2 

df 2 

<*1 

dfi 1 0 : 10:2 

dib_ 


Qi , ui 

-0i- 

1 dw~ 

+ 





_(h+t)G c 1 Ri 

«i d£i_ 





8 Q 2 

0,2 , «2 

Q 

1 dvf 

+ 





Ro 

P2 

“2 d£ 2 _ 






)> 

)> 


8<T{ m OTmH — \d\a 2 d^id^ 2 — 0 

L tic «l«2j ) 


(53) 


As all nine variations in equation (53) are independent of 
each other, it follows that the contents of all nine brackets in 
equation (53) must vanish separately. Thus the following 
nine stress-strain relations are obtained for the sandwich 
shell, indicating with an asterisk those which appear in final 
form, 


iv„r t 

, (^+f)(7*n 

iv 22 r (A+*)C*-| 



h 12 E C R{+ 

(7* 12£ c R,rJ 


J_ dwi 

u 2 d«i | 

w . (A+<)g 

(54) 

«i d^i 

0 : 10:2 d^ 2 1 

Ri 1 12 E c Ri 

N 22 n , 

(h+t)C*l 

N n r (a+oc r * 1 


c* L H 

h 12E c R 2 2 J 

<7*L" 127? c RiR 2 J 


1 d« 2 , 

Ui d a 2 | 

w (h+t)q 

(55) 

02 dfe + 

0 : 10:2 d^i 1 

R 2 1 12i? c R 2 

2 ( 1 + »') v — 011 
C* ^ 1 2 02 

A MiS A A') 

df 2 W + «i dfi \a 2 / 

(56)* 

M u — vJ \£ 22 1 

j 02 dai k 1 

(57) 

D* 

«i d^i 

aia2 CL\Ct2 (A + 0*®1 

■M-22 — vM-\\ 1 C)/?2 

j /3i ^ 1 

(58) 

D* 

0:2 dfo 

aio: 2 d^i aia 2 {h + t)R 2 


Q-2 _ R . 1 dw Ua 

(h+t)G c 2 «2 d£ 2 R 2 

G fro 7.' 1 

E c a I a 2 (h + 1 ) 


( 61 )* 

(62) 


It may be verified that the meaning of the quantities u u u 2 , 
w, 0 U and 02 is as has been indicated by comparing equations 
(54) to (61) with the corresponding equations of reference 2 
for the homogeneous shell with E c —G c = <*> . 

The system of equations (54) to (62) may be brought into 
a slightly more concise form as follows : Define the quantities 
Xi, X 2 , and X ]2 by 

, ^ (h 4- t)t E f 

Al ~2 B0 2 E c 


_1 ( h+f)t E r 
R 2 * E c 

\ 1 {h+t)t Ef 

Al2 2 R,R 2 E c J 


(63) 


and eliminate k from equations (57) and (58) by means of 
equation (62) and the equilibrium equation (34). Retain 
equations (56) and (59) to (61) in their original form and 
write for equations (54) and (55) 


Mli= ?i & (£i)+* A (ii) 

D* a 2 d£ 2 \«i/ d£i \ a 2/ 


(59)* 


( 


1 + Nn — 



N 22 =C* 


/J_ bui 

\«i 


u 2 dai . w\ . 
oil “2 d£i'Rl/‘ 


Ql a | 1 dw 

C h + t)G c Pl Aid£i Ri 


( 60 )* 


(h+t)C* 

12 E e Ri q 


(64)* 
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(i+ix,)x,,-(,.A x ,,) Nl , =0 .(Lgg 


Ul &«2 . w\ 

aia 2 d|i RJ ' 


(h + t)C* 

1 2 E c R 2 2 


(65)* 


Equations (57) and (58) become, since according to equation 
(47) D*/(h+t)E c = {mh±t)E,/E e> : • 

J_ dji . P 2 daA , H* 

«1 b|l aiQ!2 b% 2 )' E c Rx S 

(66)* 


(l + Xi)Mii — (r — \u)M 2 2= D* ^ 


(1 4'X2)-M 2 2 — (r — \u)Mu = D* 


/ 1 502 , ffl bffi 2 \ £>* 
\a 2 ^£2 ‘ aia2 d|i /' E C R 2 


s 


(67)* 


With these last transformations there is obtained a system 
of equations which is formally equivalent to the correspond- 
ing system of equations for the homogeneous shell. The 5 
equilibrium equations (25) to (28) and (31) and the 8 stress- 
strain relations (56), (59), (60), (61), and (64) to (67) are 
used for the determination of 13 quantities: Five stress 
resultants iV n , N 22 , N n , Q u and Q 2 \ three stress couples 
Mn, M 2 a, and M ]2 ; and five displacements and changes of 
slope u u u 2 , w, fa, and /3 2 . The quantity <r im which occurs 
in the sixth equilibrium equation (equation (34)) may be 
determined directly, once the shell bending and stretching 
problem has been solved. 

It is seen that the effect of transverse shear deformation 
enters equations (60) and (61) oidy and that, when G c — co , 
these equations give the values of the known theory of 
homogeneous shells without transverse shear deformation 
(references 2, 3, and 4). 

The effect of transverse normal stress deformation enters 
equations (64) to (67) only. It is seen that it is, in part, 
responsible for the occurrence of apparent stiffness factors 
(7*/(l + const. X) and Z>*/(1 -J- const. X). Thus, according to 
equation (63) , the effect of finite E c is to make the shell more 
flexible in bending and stretching than it would be with 
E c — oo . This effect, however, is present only in curved 
structures and not in plates and straight beams, as the 
quantities X have one or both of the radii of curvature in the 
denominator. A further effect of finite E c is occurrence of 
the external load terms q and s in the stress-strain relations. 
Both these effects represent, roughly speaking, what happens 
to the shape of an element of the composite shell if the length 
of the core fibers in transverse direction is changed, without 
any stretching or compressing of the face-membrane ele- 
ments. 

Having derived the general system of equations for the 
small bending and stretching of sandwich-type shells, it 
remains to apply these equations to specific problems which 
may be of interest and to determine the quantitative effect 
of the terms which are characteristic of the sandwich-type 
shell. Some of this work is done in part II of the present 
report, which follows. 

It may be stated once more that for these specific applica- 
tions the five equilibrium equations (25) to (28) and (31) and 
the eight stress-strain relations (56), (59), (60), (61), and 
(64) to (67) are used. 

I 


II— APPLICATIONS OF GENERAL THEORY 

FLAT PLATES 

The problem of the flat plate is considered first in order 
to show that the results of reference 1 are contained in the 
present results and in order to solve some problems in the 
theory of plates which have not been solved in reference 1. 

Rectangular plates. — Using notation which is customary 
in plate theory there is set 


!i = z 


a\ = a 2 — 1 

R l — R 2 ~ oo 


Ui = U 

u 2 =v 

Pl = Px 

Pi = Py 


Nn=N x 

N i2 = N X y 

n 22 =n v 

Qi = Qx 

” (68) 

II 

M n =M x 

M 12— M X y 

ML 22 ' == JMy 


V\=Vr 

P*=Pv 

m,i = m x 

m 2 — m„ „ 



The equilibrium equations (25) to (28) and (31) become 


&V* . bN xl 
by 


-\~Px — 0 


t£+!£+«-o 


(69) 


bM x dM Xi 
d.r by 


‘ Qx~ f" — 6 


bMzv | d ^-a+m„=0 


(70) 


cte by 

The stress-strain relations (56), (59), (60), (61), and (64) 
to (67) become 

bu 


N x — vNy=C* 
Ny— vNx-G* 


bx 

bv 

by 


2< 1 + ,)X.,=C*(|f+i§) 

&=(A + f)f7 e (ftr + lf) 
Q„=(h+t)G c (/3y + ^') 


(71) 


(72) 


M x - vM v = D* 


bx 

My-vM x =D* ^ 
by 

2(1 + (|f+ 5 &) J 


(73) 


As in the small-deflection theory of homogeneous plates, 
the equations for stretching (equations (69) and (71)) are 
independent of the remaining equations for transverse bend- 
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ing. Equations (69) and (71) for the stretching are not 
affected by the elastic properties of the core layer. 

Equations (70), (72), and (73) have been treated in refer- 
ence 1 by means of a stress function \}/, which, together with 
the deflection w, was taken as one of two basic variables. 
In what follows an alternate treatment is given, in which the 
problem is reduced to three simultaneous equations for the 
quantities ft,, /3„, and w. On the basis of these three simul- 
taneous equations a problem not considered in reference 1 
is treated, namely, the bending of a rectangular plate which 
is simply supported on all four edges. This same problem 
has also been solved by Donnell by a method which differs 
from the one employed here. (See reference 5 where the 
case of the homogeneous plate is considered.) 

To reduce equations (70), (72), and (73) to three simul- 
taneous equations for p x , /3„, and w, first a quantity co is 
defined by 




(74) 


Introducing equation (72) into the first of equations (70), in 
view of equation (74), there is obtained 

a+X ,2 w=—<i/(Ji+t)Gc (75) 

where V 2 — d 2 /dx 2 -\-d 2 fdy 2 . 

Next, Q x , M x , and M xy are taken from equations (72) and (73) 
and the result is substituted in the second of equations (70). 
This gives, after slight transformations, 


~ v ^ x -2(h+t)G c p x +^^-2(h+t)G c wj+m x ==0 

(76) 


In an analogous manner the following further equation is 
obtained : 


^ V 2 0- 2(h+l)G c 0 v +^ 2(A + f) &«?]+ m,= 0 

(77) 

In order to solve equations (75) to (77) two equations are. 
next obtained involving w and co only. Differentiating 
equation (76) with respect to x and equation (77) with respect 
to y and adding the two resultant equations, in view of equa- 
tion (74), give 


The following procedure may now be carried out: (a) 
Solve equation (78) for co, (b) with this value of co solve equa- 
tion (75) for w, (c) substitute co and w in equations (76) and 
(77) and solve for (i x and p„, and (d) eliminate extraneous 
terms in (3 X and /J„ by considering equation (74). 

Before deriving the solution of a problem along these lines, 
the explicit differential equation for w which follows by 
combining equations (75) and (78) may be given 


Note that the effect of transverse shear occurs on the right 
side of the equation only. In order to compare the magni- 
tude of the q terms on the right of equation (79), assume that 
relevant changes of q occur over distances of order l (where 
l may or may not be a representative diameter of the plate) . 
Then, as order-of-magnitude relations, there result 



V 9 -q 

( h-\-t)G c 


0 



j 


(80) 


From equations (80), it follows that transverse shear ceases to 
be a secondary effect as soon as l is of order pht V E f jG c or of 
smaller order. 

Bending of rectangular plate with simply supported 
edges. — The edges of the plate are assumed to be at x=0,a 
and y=0,b and along these edges moments and deflections 
are assumed to vanish. Further, 

m x =m y = 0 

CO CD 

g=Z X) imn sin \ m x sin y n y 

m = 1 7i’= 1 

where 

X m = rmr/a\ 

jU TC = 7l7r/& ) 

From equation (78), it follows that 

C0=-L SE — a sin b m x sin Auy + co* (83) 

\ m Tfl« 


(81) 


(82) 


on* 

~— 2 V 2 w-2(h + t)G c (a+ V 2 w) 

1 — V 


dm x bm y 
c)x c )y 


0 


where w h is a harmonic function. Putting equation (83) into 
equation (75), 


and, making use of equation (75), 


V 2 co=- 


D 



C)TOy \~] 

*v ) J 


V 2 w 




HD 


-t)G c _ 81 


(78) 


_Xto 2 +m* 2 (h-\-t)G c 
which is integrated to 


sm \ m x sin n n y — oj h 




Qm n 

(bJ+Vn 2 ) 2 


[‘ 


P(Xct 2 +M» 2 ) ~] 

(h+t)G e J 


sin \ m x sin fi n y+w h 


(84) 


where w h is the general solution of V 2 w n — — co h . It is to be 
expected and may be shown explicitly that for the plate 
which is simply supported all around w h — co A =0 and, as in 


the Navier solution for the plate without transverse shear 
deformation, the particular integral is the complete solution 
of the problem. 
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Equation (84) may be rewritten in the more explicit form 


w- 


' D 7T 4 


2S 


»[i 


E f 


2 (1 -v 2 )G c 


^±^(m 2 +?i 2 


£!Y| 

ft*/!: 


[m 2 +-?i 2 (a 2 /6 2 )] 2 


mir . 

sm x sin — j~ y 


11 


(85) 


When G c = c° , equation (85) reduces to Navier’s solution 
w N . Equation (85) is more readily interpreted by means of 
the ratio w/% of deflection with and without transverse 


shear deformation. On the basis of equation (85), there may 
be obtained the following equation (86), ' 


DV 2 w n 

w Wn (h+t)G~ c _ (h+t)tE f V 2 w n 

w N w N 2(1 — V*)G C w N 


= 1 + 


7 r 2 Ef (h-{-t)t 

"2" (1 — v 2 )G c a 2 




sin mirx/a sin mry/b 
m 2 +n 2 (a 2 /fr *) 




q mn sin rmrx/a sin rnvyjb 
[m 2 +n 2 (a 2 /b T )] 2 


( 86 ) 


Setting -ir 2 /2 (1 — v 2 ) — 5A and (E,jG c )t{h J r t)/a 2 =(S, equation 
(86) takes on a form which contains as a special case Donnell’s 
result of equation (18) of reference 5. 

For the case of a uniform load intensity q = Constant and 
for the center of the plate (x—a/ 2, y=b/2) equation (86) 
becomes 

sin rmr/2 sin rnr/2 

w _ , tv 2 E f (h + t)t mn[m 2 +n 2 (a 2 /b r )] . . 

w N '2 (1 — v^Gc a, 2 sin tott/ 2 sin nir/2 ' 

mn[m 2 -\-n 2 (a 2 /b 2 )] 2 

The ratio of the series is 1.98 when a/b= 1, and the ratio of 
the series is 1.11 when a/b = 1/2. 

For the case of a concentrated load at the center of the 
plate the deflection ratio at the point of load application 
assumes the form 

(sin mw/2 sin m r/2) 2 

w 7T 2 Ef m 2 +n 2 (a 2 /b 2 ) 

w N 2 (1 — v^Gc a 2 -sr-'-sr-' (sin wiir/2 sin mr/2) 2 ^ 

2 - j2 -' [m 2 +n 2 (a 2 /& 2 )] 2 


Now it is easily shown that the numerator series in equation 
(88) does not converge and consequently w/w N = in this 
case. A more detailed consideration shows that in any plate 
theory which takes transverse shear deformation into 
account the deflection under the point of application of a 
concentrated load must become infinite in contrast with what 
happens when transverse shear deformation is not taken into 
account. This difference, of course, vanishes as soon as the 
load intensity becomes finite, and then the theory with trans- 
verse shear deformation taken ’into account is more accurate 
than the theory which does not take into account this effect. 

For the sake of numerical illustration take again the 
square plate (ajb= 1) with uniform load distribution. 
According to equation (87), the deflection at the center is 
increased because of transverse shear by the factor 


w 

Wn 


= 1 + 9.7 


Ef 

G c a 2 


(87a) 


Take A=1.0 inch, f=0.1 inch, a=10 inches, E//(? c =200, and 
r=l/3. Then, according to equation (87a), w/w. v = 1+2.3, so 


that in this case the deflection with transverse shear is more 
than three times the deflection when shear deformation in 
the core is neglected. 

Returning now to equation (84) for w and equation (83) 
for to and substituting these two equations in equations (76) 
and (77) in order to determine the changes of slope dr and fi v , 
after slight transformations there results 


dr=— 7) g r? x cos ^ sin 

p 9 =— Vj ES q rZ >ln 2\2 sin C0s VnV 

U \A-rn ~T~ A * n ) 


(89) 


Equations (89) are remarkable for the reason that they are 
not affected by transverse shear deformability. According to 
equations (73), the same is then true of the bending and 
twisting couples M z , M v , and M xy . It is not easy to see 
why, in this statically indeterminate problem, the magnitude 
of the internal forces does not depend on the elastic proper- 
ties of the core. The analysis, however, shows that the 
distributions of M x , M y , and M xv , and therewith, of Q x and 
Q y , remain the same as those obtained under the assumption 
that G c = co . In this connection the following remark may 
be made. 

Evidently the following three boundary conditions, 
w—M z =p v ~ 0 along the edges x=0,a, have been satisfied. 
In order that the last of these three conditions be satisfied 
there are necessarily nonvanishing edge values of the twisting 
couples M xv . The same is true in the theory without trans- 
verse shear deformation, where, however, no alternative 
possibility exists; as in that theory only the boundary condi- 
tions w=M x = 0 are relevant. For the present system of 
equations three boundary conditions must be formulated for 
every plate edge. Thus, it is possible although mathemat- 
ically complicated to solve the problem of the rectangular 
simply supported plate with the edge condition P v =0 
replaced by the condition M xv — 0. In that case, which will 
not be pursued here, there evidently will be a distribution of 
internal stresses which is modified by the effect of transverse 
shear deformation. 


888244—50 3 
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Cylindrical bending of plates. — As a further relatively 
simple example of application of equations (70), (72), and 
(73) problems are considered for which 


d( )/dy = 0 

d( )/c)x~d( )jdx = ( )' 
'hd X y= Qy~ 7Yly= (3 l 0 
M u = vM x 


(90) 


and where consequently the problem reduces to the following 
system of equations : 


Qx -\-<i= o 
M x ' — Qx+m x = 0 
Qx = (h + 1 ) G c (p x + w') 
(l-v*)M x =D*p x ' 


(91) 


Note that in order to obtain the problem of the sandwich 
beam from equations (90) and (91) the only changes which 
are necessary amount to setting r=0 in equations (91). 

To set into evidence the effect of finite values of G c in 
equation (91), the following system of equations is deduced 
from equation (91) : 


DW’-q+m,' 

(92) 

«•= D ”>" (*+i )ft 

(93) 

Q x - Dw { j l j r{)G +m x 

(94) 

^~ w ' + (h+t)Gc 

(95) 


Solutions to the following problems are listed: 

(1) Simply supported plate oi span l carrying a load 
q- --r/ 0 cos ttz/L Boundary conditions: w(±l/2)~M x (±l/2) 
= 0 . 


q a F, , r 2 E f cos -irx/1 

W ~D L + ~2 (l — v 2 )Gc J (tt/0 4 " 


(96) 


The factor in brackets may again be written in the form 
1 +5.4)3, with p=(E f /G c )/[(h+t)t/l 2 ], using the notation 
suggested in reference 5. As the problem is statically 
determinate as far as moment and force are concerned there 
is no modification of M* and Q x due to the finite value of G c . 

(2) Simply supported plate oi span l carrying a uniform 
load g=go ■ 

From this there is obtained for the center deflection, 


w(0) = 


5 q 0 l 4 
'384 D 




1 + 


24 


E 


5 (1- 


(A-f-t)t~ 
v*)G c G _ 


(98) 


It is seen that the correction factor for the center deflection 
is almost the same as that for the cosine load curve (equation 
(96)), the only difference being a change of the factor 
ti 2 /2=4.93 into 24/5=4.80; that is, a reduction of the shear 
correction factor by at most 3 percent is present. Note 
that according to equation (87a) the shear correction factor 
for the square plate of width a— l is more than twice as 
large as the shear correction factor for the plate strip of 
width l. 

(3) Built-in plate, of span l carrying a uniform load g=g 0 . 
The boundary conditions are: w(±l/2) = p x (±l/2)~0 (and 
not w' ( ± 1/2) = 0) . 

™-md ![(!)'- 0 “ 2 [ i + » lw ][ te )’- 0 ! <99) 

From this there follows for the center deflection, 

W (Q) — go - — Tl + 24 — (100) 

384D|_ + (l-* 2 )^ P J { ' 

Comparison of equations (100) and (98) shows that for 
the built-in plate the effect of transverse shear deformation 
is very much more pronounced than it is for the simply 
supported plate, a factor 24/5 in the latter case being replaced 
by a factor 24 in the former case. A somewhat similar per- 
centage increase must take place in going from equation (86) 
for the rectangular plate with all four edges simply supported 
to a formula (which has not yet been derived) for the rec- 
tangular plate with all four edges built in. As a further 
result in this problem of the built-in plate, by putting equa- 
tion (99) into equation (93), it is found that the moment 
function M x does not contain any terms depending on the 
effect of transverse shear deformation. This again is some- 
what surprising as in this case it is not possible to determine 
the moment function by statics alone. As a problem where 
the moment distribution is in fact dependent on the effect 
of transverse shear there may be mentioned the problem of 
the cylindrically bent plate with both ends built in, which 
carries a load gi=g\X instead of the load g 0 =3 o- This 
problem also may be solved by means of equations (92) to 
(95). 

Circular plates ; rotational symmetry. — As no examples of 
solutions of circular sandwich-plate problems have as yet 
been published and as it is of some interest to determine in 
which way the shear correction factors change in going from 
a problem for the plate strip to the corresponding problem 
for the circular plate, the equations for axisymmetrical trans- 
verse bending of circular plates are briefly discussed. 

Polar coordinates r,6 are introduced and notation which 
is customary in plate theory is used. As a consequence of 
equations (70), (72), and (73), the following system of equa- 
tions is obtained : 


drQr/dr-f-rg — 0 
drM T ldr—Me—rQ T +rm T =0 


( 101 ) 


Q T = (h+t)G c (p r +dw/dr) 
M r — vM e = D *d(3rldr 1 
Mo—vM r =D*B r /r ) 


( 102 ) 

(103) 
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there results 


According to equation (79), the equation for the deflection 
w will be 

V 2 q 


.11 drm T 
DV 2 V 2 w—q+ 7 z- 


2 r dr (h-\-t)G c 


(104) 


where V 2 =(l/r)d[r d{ )/dr]/dr. Having found w by means 
of equation (104), /3 r may be determined from 


0T=- 


dw 


Qr 


dw , 1 


dr ( h-\-t)G c dr 


dr 


(105) 


and therewith M r and Me are obtained from equations (103). 

In the present problem it seems to be somewhat more 
convenient to proceed as follows: Combine equations (101) 
and (103) to obtain as equation for the change of slope p r , 


r _ d 

(1 dr ft A” 

_ dr' 

\r dr ) J 


dr r 


r dr 


(106) 


Having P T , M r and Me are found from equation (103) and 
Q r , from the second of equations (101), 




Mi 


= D 0+-it) 


U-r dj £)+™’ 


(107) 


(108) 


Finally, with this value of Q r , w is found by integrating 
equation (102), 

ro j , D drfi r \ , f m r d T ,, 
w- fPrdr + ^ +()0 ^~ dr )+- {h+m (109) 

Deflection of circular plate with built-in edge. — The 

bending is now considered of a plate with transverse load 
q—q n (r/a) n and with m r = 0. First, from equation (106), 


Dp^e^+et (jV+^log^- 

a \a/ a a 


q [ntf 


(«,+ 4)(n+2) : 


er 


( 110 ) 


Attention is restricted to complete plates with no con- 
centrated load at the center, and consequently it is necessary 
to set c 2 =c 3 =0 in equation (110). This gives 


Dp r =ci - 


q n a 6 


a In -(- 4) in -(- 2) 


»+3 


(1 10a) 


Putting equation (110a) into equation (109), there results 
for the transverse deflection w 


a 12 \a / 


q n a 6 


(n+4)(n+2) 2 


GTV] 


+ 


D 


a\h + i)G, 


[^+js^en <ni) 


Taking the case of a plate of radius a with built-in edge, 
that is, with the boundary conditions 


$,(a)=w{a) = 0 


( 112 ) 


Dp r 


and 


qna 6 

'(7i+4)(7i+2) 2 


a) aj 


T-i qna* ( (r/a)^ 4 -! I {rjaf-l 

uw (n+2) 2 ( (n+4) 2 2 n+ 4 


D 


a 2 (h + t)G c 


[©'"-I 


13 


(113) 


(114) 


From equation (114) there follows for the deflection at the 
center of the plate 


D w(0)~ 


q n a 


2(m+2)(?i+4) : 


(rc.+ 4) 2 E f (h+t)tl 
_ 1+ n+ 2 (1 — v*)G c a 2 J U 


Consider the following special cases : 

(1) Uniform load distribution g„=go ■ From equation 
(115), there follows for the ratio of deflection with and 
without transverse shear deformation 


w( 0 ) 


[w(0)] , 


= 1 + 8 


Er 


(h~\~ f)t 


(1 -v 2 )G c 


(116) 


Equation (116) may be compared with equation (100) for 
the deflection of the infinite plate strip of width l with built-in 
edges. Setting l— 2a, it is seen that, while the transverse 
shear correction factor for the strip has a value 6, the corre- 
sponding factor for the circular plate is 8. This is consistent 
with the earlier comparison between the simply supported 
strip and the simply supported square plate, except that 
there the change is from 4.8 to 9.7. 

(2) Linearly increasing load distribution g=g x rja. From 
equation (115), it follows that 


w(0) 


©(0)1, 


= 1 + 


25 


E f 


(k+f)f 


3 (1 — v 2 ) G c 


(117) 


showing that the correction effect is only slightly greater 
than in the case of the uniform load distribution. 

(3) Load increasing linearly from edge to center, g=g 0 + 
gi(r/a). (</, = — g n ) . From equation (115), it follows by 
superposition that 

n g ° a4 Tl 4-8 Ef (A + © ~j 

Dw(0) 2X2X16 L 1 + 8 (l-v^Gc a 2 J 

go a* f, ,25 Ef (h+t)t l 

2X3X25 L 3 (1 -v*)G c a 2 J 

43g 0 a 4 r, ,3000 E f (h+t)f] n , cX 
L + 387 (1 — v 2 )G c a 2 J 


D w(0) 


32X150 


Comparing the factor 3000/387 = 7.76 which occurs in equa- 
tion (118) with the corresponding factors 8 and 8.33 in 
equations (116) and (117) it is seen that, in the foregoing 
three problems at least, there is little difference between the 
transverse shear stress correction factors in the case of three 
different loading conditions for the circular, clamped-edge 
plate. The fact that this agreement should not be expected 
to hold generally follows again by considering the case of a 
point load at the center of the plate, for which the shear 
correction factor would again be infinite. 



14 


REPORT 975 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


The examples of this section should be augmented by the 
solution for the circular plate of radius a, which carries a load 
distributed uniformly over a smaller circle which is concentric 
with the boundary of the plate. 

CIRCULAR RINGS 

As the simplest example of a curved sandwich structure 
there are considered in this section stresses and deformations 
of circular rings in their own plane. As was found in the 
general developments of part I of this report, in a curved 
sandwich structure there will be the effect of both transverse 
shear and normal stress deformation. (The effect of trans- 
verse shear stress deformation on homogeneous circular rings 
has been considered by Beskin in reference 6.) 

There are set for the relevant coordinates and variables 


h =ad 

a x =l 


U x = V 

0i=0 

Nu=N 

o> 

il 

M U =M 

Pi=P 

m x = m 

a( m= 

■d( )/add=( )'/a 


Eing sector acted upon by end bending moments. - : As a 

first problem on circular rings, which illustrates the effect of 
transverse normal stress deformation, there is taken this 
basic case for which, as is known, there must be the same 
stress distribution at all sections 8 = Constant of the ring. 
According to equations (120) and (121), 


N=Q = 0 

M=M 0 f 

c T tn = —M 0 /(h+t)a , 


(126) 


Equations (122) to (124) become 


v' ~\-w= 0 
j3 + (w'—v)/a = 0 


(127) 


(1 + X)M 0 = (128) 


The significant result of this consideration is contained in 
equation (128), which may be written in the alternate form 


Z>* 

M=M 0 =— 

a 


1+A~ 


1 +-. 


D* 

Ef 

a 2 E c 


w"-\-w 


(129) 


Xl =A =\[(h+t)tla 2 ]{E r IE c ) 


The equilibrium equations (25) to (28), (31), and (34) 
reduce to the following equations: 


N' + Q + ap^0^ 

Q' —N+aq = 0; 


( 120 ) 


M'—aQ+am = 0 


( 121 ) 


crtm = s—[M/(h-\-t)a] ) 

The stress-strain relations (56), (59), (60), (61), and (64) 
to (67) reduce to the following equations: 

( 1 +I x ) W= S < 7 *[''+“’+ ! I5#] (i22) 

Q = (h+t)GJf}+±(w'-v) ] (123) 

(1+A )M=i£>*(V+^) (124) 

The load terms p, q, m, and s are given, according to 

equations (19) to (22), by 


P = 

O+^'M 

h+t\ 


2 = 

0 +'£>*+( 

h-\rt\ 

1 2„> 

►- 

m= 

W+&0 

* i 1 \v> 

■] 

s = 


-0 



Thus, in this case of pure bending the transverse flexibility 
of the core is responsible for a reduction of the bending stiff- 
ness factor D*=(}{)t(h J ct) 2 Ef which is obtained exactly when 
E c = 0 and practically when E c is of the same order of magni- 
tude as E f . Equation (129) shows that the reduction of D* 
is significant whenever E c is so small that the ratio EJE f is 
of the same order of magnitude as the ratio ( h-\-t)t/a 2 . 

As a numerical example take the following values: A=0.9 
inch, f=0.05 inch, a=20 inches, and E f /E c =1000, for which 

l»i^|_‘ X M^ 06 xlo oo = 0 . 0 595 


indicating a reduction in bending stiffness of about 6 percent. 
Changing a from 20 to 10 inches changes the effect from 6 to 
24 percent. Changing Ef/E c from 1000 to 2000 increases the 
effect from 6 to 12 percent. Altogether it may be said that 
this effect is of noticeable magnitude for some geometrically 
reasonable structures when the modulus ratio E f /E c is of the 
order 1000 or more. Assuming aluminum face layers with 
E f = 10 7 psi, this means that E c « 1 0 4 psi, which is well within 
the range of some present-day core-layer materials. 

Comparing equation (129) with the earlier formulas for 
the effect of transverse shear stress deformation, for instance 
with equation (116) in which a represents the plate radius, 
and observing that G C ~(Y 2 )E C , it is seen that the correction 
terms are of the same form, the difference being an appreci- 
ably larger numerical factor in the expression representing 
the shear effect. 

Closed circular ring acted upon by uniform radial load. — 

Having rotational symmetry, d/dd = 0 and r = 8 = 0 . Also set 
p=m=0. The remaining equations permit the determination 
of the stresses in the face and core layers in a way which 
depends on the extent to which the load is applied to the 
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outer (upper) and inner (lower) face membranes. Equation 
(12) becomes 

N=aq (130) 

From equations (121), it follows that 

<r{ m =s—[M/(h+i)a] (131) 


The stress-strain relations' (122) to (124) give 


w 1 A , X\ 

e-^+s) 02 ' 


( h+t)q aq 
l2E c a~~C* 


r, , * (h+t)C * “I aq 

L ‘‘"a 1 2 E c a 2 J C* 


and 


M— 


D* s 
1 + A a E c 


(132) 

(133) 


A closed circular ring subjected to a uniform radial load 
distribution q is stressed not only by a uniform axial force 
N=aq, as would be expected, but in addition is stressed by 
a uniform bending moment M, the magnitude of which is 
given by equation (133). The explanation of this result is 
that for a ring with relatively soft core the circumferential 
stress distribution depends on the extent to which the 
external radial load is applied to the inner and outer forces, 
respectively. Roughly speaking, for a sufficiently flexible 
core layer the load q u goes predominantly into the outer 
face layer, while the load </, goes predominantly into the 
inner face layer. 

According to equations (46), the stresses in the two face 
layers are given by 


N -=* N +hTt M 

N ‘=2 N -kTt M 


(134) 


According to equations (130) and (133) and in view of the 
definitions of D* and A, this may be written 



(135) 


Combining next equations (131) and (133), for the trans- 
verse normal stress in the core layer, the following expression 
is obtained: 


m 


s 

1+A 


(136) 


For a specific example assume that the radial load is 
applied entirely to the inner face of the ring so that q u — 0 
and, according to equation (125), 




h+t\ 

2 a J 


<h 


With q and s given by equations (137), equations (135) and 
(136) become 


■ (l h+t " 

\ m 

V 2a , 

f l+X 

► 

/, A+t N 

\ (1 d-2X)ag r ! 

V 2a y 

> l+X 

- Ul- 

A+f\ qi 

2 \ 

2 a J l+X 


(138) 


(139) 


It is seen that the flexibility of the core layer increases the 
circumferential stress in the loaded face layer in the ratio 
(1 -)— 2X) /(I -j— X) and decreases it in the unloaded face layer 
in the ratio 1/(1 +X), where X is defined by equation (119), 
compared with the equal values of these stresses when 
E c = co . 

Considering once more the numerical data under the sec- 
tion entitled “Ring sector acted upon by end bending 
moments,” it is found, for instance, that the stress in the 
inner face layer may be about 6 or 12, or 24 percent higher 
than the corresponding stress calculated without taking into 
account the transverse flexibility of the core layer. 

Ring sector acted upon by radial loads q u and q h uniform 
in circumferential direction and with vanishing resultant q . — 
Again it is assumed that d{ )/dd= 0, m=p= 0 and now in 
addition that g=0, so that, according to equation (125), 
the only nonvanishing load term is s. Further, it is assumed 
that the ends 8— ± a of the ring sector are free of stress, 
that is, Af(±a) = f?(±aO=M(±a:) = 0. The ordinary theory 
of circular rings would then indicate the absence of deforma- 
tions in the entire ring. In the present case there is found 
a type of deformation peculiar to the sandwich ring, which 
may perhaps be compared to the action of a Bourdon gage. 


Solving first equations (120) and (121) and satisfying the 
end conditions of the ring sector, 


N=Q=M = 0 

Cfm = S 


(140) 


The stress-strain relations (122) to (124) are then 
v '+w= 0 

— v — 0 ► 

0'— — s/E c . 


(141) 


Assuming s independent of 8, from equation (141) there is 
obtained by integration, with constants of integration A u 
A 2 , and A 3 , 


0 = — 


v = —a~-d+Aia + A 2 cos 0+A 3 sin 8 
tie 


( 142 ) 


w=a^r+A 2 sin 8— A 3 cos 8 


(137) 
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As a specific example consider a complete ring, slit ted 
radially at the section 6=ir, so that a—ir. Prescribe fur- 
thermore the symmetry conditions /3(0) =»(0)=w(0)=0. 
Under these conditions there is obtained from equation (142) 


Ej3=—s8 

E c v=—as(9 — sin 6) *- 
E c w=as{ 1— -cosfl) „ 


(143) 


From equations (143), it follows that the radial slit, which 
is of zero width before the loads q u and q t are applied, opens 
under the action of the loads to a width given by 

»(— ») — v(%) = 2ira -^-=2ira ^1 ^ (144) 

For a numerical example take a= 10 inches, h — 1 inch, t= 0.05 
inch, E c = 10,000 psi, and g«=20 psi, and obtain 

v( — x) — -v (x) = 0.132 inch (145) 


± («'-»)=& sin 

— n ■ n a 2 r . I 1 (h+tyt ( Ef , -E/M n rns 

-Qosin ^p*[_l+ 2 a 2 \e + gjJ ( 15 °) 


Simultaneous solution of equation (150) and the first of 
equations (148) for v and w gives as general expressions for 
v and w, 


v — A8 cos d-\-Ai sin 0+ A 2 cos d 
w=Ad sin 8 — {Ai^rB) cos 8AA Z sin 8 


(151) 


where A i and A 2 are arbitrary constants of integration and 
A and B are found to be 


a-m. n .a ( E r i e a i p* (i 4>yi 

^~2Z9*L + 2 a 2 \E C ' G c )'~ a 2 C* \ ~^3/J 

l {hAt)t(E f E f \ D* / \\- 1 

^ _ 2L>*L + 2 a 2 \Ej GJ o*C* \ 3/! 


r (152) 


The foregoing three examples of ring analysis have been 
discussed in some detail, because they illustrate relatively 
simply the effect of transverse normal stress deformation in 
the theory of curved sandwich structures, without involving 
at the same time the effect of transverse shear stress 
deformation. 

Bending of semicircular ring by end shear forces. — A 
problem is now considered in which both the values of E c 
and G c affect the result of the analysis. In the equilibrium 
equations (120) and (121) all external load terms are set 
equal to zero and then, by integration and from the boundary 
conditions, that is, from 


iv(±f)=M(±f)=° 

Q -0= ± Qo 


(146) 


The following expressions for N, M, and Q are obtained : 


Q=Q 0 sin8 

N=Q 0 cos 8 ► 

M——aQ 0 cos 6 „ 

The stress-strain relations (122) to (124) become 
(l+A/3)Qo cos 8=(C*fa)(v' + w) 

Qo sin 8=(h+t)G c [l3 J r(w' — v)/a] ► 
— (14-X)Q 0 ® cos 8=(D*/a)&' 
Integration of the last of equations (148) gives 
Z?*/3 = — a 2 (l -f X) Q 0 sin 8 


(147) 


(148) 


(149) 


As further conditions, it is prescribed that »(0) —v(ir/ 2) =0, 
which makes A 2 — Ai = 0 in equation (151). There remains 


v = Ad cos 8 

w= Ad sin 8—B cos 8 


(153) 


Of particular interest are the values of w(ir/2) and w(0), 
the first of these giving the radial deflection of the point of 
load application, the second giving the change of radius at 
right angles to the applied load. It is found that 


(154) 


, n s_ 1 Qo® 3 

w{0 ) — B 2 j-y* 




1 + 


1 Ef . Ej 


\Ec G c ) a 2 G*\ i+ 3;J 

(155) 


Equations (154) and (155) contain the interesting result 
that, for this problem, transverse shear and transverse normal 
stress affect the outcome formally in nearly the same way. 
If the generally unimportant terms with D*/a 2 C* are omitted, 
which amounts to the usual assumption of circumferential 
inextensibility of the ring, then the effects of finite E c and 
G c occur in exactly the same way. 

For a numerical example take h= 0.9 inch, t=0.05 inch, 
a— 20 inches, E f /E c — 1000, and E f /G c = 2000. This gives 

D* _1 (h+t) 2 _ 1 

a 2 C* 4 a 2 1770 

1 (h+t)t E f 2 

2 a 2 G c 16.8 


1 (h + t)t_ 1 

2 a 2 ^16,800 

1 (h+t)t Ef 1 

2 a 2 E~ 16.8 


The factors in brackets in equations (154) and (155) become 


14 

14 


1 

16.8 

1 

16.8 


*“16.8 

f _2_ 
r 16.8 


1 

1770 

1 

1770 


Q+sxio)- 1 18 

( 1+ 3Xl6l)” 118 


where a constant of integration has been eliminated by means 
of the symmetry condition /3(0) = 0. Substituting equation 
(149) in the second of equations (148), 


and 
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Thus, in the present example the flexibility of the core is 
responsible for an 18-percent increase of deflection-load ratio, 
and of this 12 percent is due to transverse shearing and 6 
percent to transverse normal stress. Compared with these two 
effects the effect of circumferential extensibility of the com- 
posite ring is seen to be negligible. As a further numerical 
illustration, it is noted that reducing the ring radius a from 
20 to 10 inches, with all other data unchanged, changes the 
18-percent correction to a 72-percent correction. 

Bending of complete circular ring under action of two 
concentrated radial forces at 6— ± 7r/2. — The solution of 
this problem may be obtained by superposition of the solu- 
tions for the semicircular ring under the action of end shear 
forces Q 0 (equations (146) to (155)) and under the action of 
end bending moments M 0 (equations (126) to (129)). 

The first step consists in determining M 0 in terms of Q 0 
such that the sum of the /3's from equations (129) and (149) 
assumes the value zero for 0=x/2; that is, the value of the 
superimposed bending moment at 0=ir/2 must make the 
tangent to the deflected ring at this point horizontal. Com- 
bining equations (129) and (149) in this manner, there is 
obtained 

2 jy* o 0 

or 

M u =(2/ir)aQ 0 (156) 


It may be noted that equation (156) is a further case of a 
statically indeterminate problem where transverse shear and 
normal stress flexibility do not affect the internal force and 
moment distribution but affect only the state of deforma- 
tion of the structure. 

Further, the radial deflections w(-ir/2) and w(0) due to the 
action of M 0 are calculated, in order to combine them with 
equations (152) and (153). Integrating equations (129) and 
(127) with the boundary conditions ®(0) =v(w/2) =0, there 
is obtained for the displacements due to M 0 , 


D*w— — (1 +X)Af 0 a 2 ^l — ^ cos 6 ^ 

D*v =(l + X)M 0 n z ^-|sin 6^ 
and, in particular, 

D*w (0) =(l+A)M 0 a 2 (|-l) 

> 

T»*w(|)=- (1 + X)M°a 2 


(157) 


(158) 


Combining equations (158) with equations (154) and (155) 
and taking M 0 from equation (156), there follows for the 
resultant displacements 


w 



Qod? 

D* 


w{ 0) 


Q 0 a 3 
D* 


{(!“?) (!+ X )+j[ X « 

}0 _ 1) (1+X)+ I[ X ° 


D* 


a 2 C* 




X (159) 


where \ G 


^ -vr has been put as a further abbreviation. 
2a 2 G c 


Equations (159) may be written in the alternate form 
w(j)=0.149^- 3 |l + X + 5.29[x G +^ 5 (l-l-|)]| (160) 

w(0) = -0.137^|l+X + 3.65[x o -^(l+!)]| (161) 

When X=X G =0 and when the composite ring is assumed 
axially inextensible, which amounts to putting D*/a 2 C *= 0 in 
equations (160) and (161), then equations (160) and (161) 
reduce to well-known results of circular-ring analysis. 

Comparing equations (160) and (161) for the closed circular 
ring with equations (154) and (155) for the open semicircular 
ring, it is noteworthy that for the semicircular ring X and X G 
occur with equal weight, while for the closed circular ring the 
influence of X G is considerably greater than the influence 
of X. Thus, for the closed circular ring the effect of 
transverse shear deformation is much more important than 
the effect of transverse normal stress deformation, while 
for the open semicircular ring both effects occur in a much 
more nearly equally important way. 

For a numerical example of the use of equations (160) and 
(161) take again the values for the numerical example given 
in the section entitled “Bending of semicircular ring by end 
shear forces.” This gives for the expressions in braces 


11 1 J 

, 2X5.29 

5.29 

1 1 16.8 

h 16.8 H 

1770 

1 | 1 i 

, 2X3.65 

3.65 

1 + 16.8 H 

1 16.8 “ i 

1770 


Thus, while the effect of transverse stress deformation for the 
open circular ring amounted to 18 percen , the corresponding 
corrections for the closed ring are 69 and 50 percent, 
respectively. 

The next step in the analysis ol sandwich-type circular 
rings would be the general solution of the system of equations 
(120) to (124) for arbitrary load distributions. This, evident- 
ly, is possible and further specific examples of interest might 
be analyzed on the basis of the general solution. 

CIRCULAR CYLINDRICAL SHELLS 


In this section the general system of equations of part I 
of this report is restricted to the equations of the theory of 
circular cylindrical shells. The treatment of sandwich-type 
shells of this kind is shown to be not appreciably more 
difficult than the analysis without the effect of transverse 
shear and normal stress. 

As specific examples, some problems of rotationally sym- 
metric deformations are treated. In particular the influence 
coefficients are obtained for a semi-infinite shell acted upon 
by bending moments and transverse forces at one end of the 
semi-infinite shell. With these influence coefficients an 
explicit solution is obtained for the problem of the infinite 
circular cylindrical shell acted upon by a pressure band of 
zero width. 
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In the general equations of the problem there are set for 
the relevant coordinates and variables, 


=ad 

&=x 

Bi=a 

H 2 = CO 

Nu=Ne 

N 2 2 = N x 

Qi~ Qe 

Qi = Qx 

Mu=Me 

M 2 2 — M x 

Ux = V 

U 2 ~U 

/3i==/3e 

II 

cq 

mi=m e 

m 2 =m x 

Pi—pe 

Pl=Px 


ati — ot2~ 1 



(162) 


iV"i 2 — N 2 l — N x 0 -1/(2 — l/v ( — -l/r; 


The equilibrium differential equations (25) to (28), (31), 
and (34) become 


d/Vs , 1 dN xt 
dx 

dN xS 


a c )8 


-\-Px— 0 


1 dN» Q, 
a 58 a 




dx 

dx ^ft d0 a " 1_a u 


(163) 


dMj , 1 d Mr, 


dx ft d0 
dMrf 1 dM« 


dx a de 




“f" Wlx — 0 
■ Qe-\- Me=0 
Me 


0 h-\-t)a 


(164) 


(165) 


The stress-strain relations (56), (59), (60), (61), and (64) to 
(67) become, with X 2 = Xj 2 = 0 , A, = ( 1 /2) [(A + t)t/a*] ( E f /E c ) = X , 


H) 


N x - 


Ne 

vNe=C* ( 


..AT —n* f 1 £>» , w , (h+t)q 

vNx -° lade + a + l2^E c ] 


d u\ 
dx/ 


2(1 + W--0*(g+|g) 
e,=(*+ea(ft+|||-|) 

Q.-(h+()G, («.+§) 


1 (166) 


(167) 


._f_A 

aE c ) 


2(1 + -».-»* (H'+i^) 


( 168 ) 


When G C =E C = <» (and therewith X=0) equations (163), 
(164), (166), (167), and (168) reduce to the known system of 
equations in which deformations due to transverse stresses 
are neglected. The solution of the present system of equa- 
tions is not essentially more difficult than the solution of the 
system with G C =E C = =° . In particular also here there may 
be obtained a trigonometric double-series solution, as a 
generalization of Navier’s solution for the flat plate (references 
7 and 8). 

For this trigonometric double-series solution there is set, 


q — X X imn sin md sin nx/l 
Po = X ! X Pemn cos md sin nx/l 
Pz=X X Pxmn sin md cos nx/l 
m x =X X m xmn sin md cos nx/l 
m<)=X X mamn cos md sin nx/l 
<=2S*« sin md sin nx/l 


(169) 


(170) 


(171) 


w= X X. w mn sin md sin nx/l ' 

« = X X cos md sin nx/l 
w = X X u mn sin md cos nx/l >- 
ffi— X X) sin md cos nx/l 
fie = X X P»mn cos md sin nx/l ^ 

Q Z =XX Qxmn sin md cos nx/l 
ft=XX Qomn COS md sin nx/l 
(N x , iV*) = XX /V« m „) sin md sin nx/l 

A(i=XX iVzemre cos md cos nx/l 
(M x , Me) = x X (M zmn ,Me mn ) sin me sin nx/l 
M r e = X X Mximv. cos md cos nx/l 


When equations (169) to (171) are substituted in equations 
(163) to (168) there remains for every value of m and n a 
system of 13 simultaneous equations for the 13 Fourier 
coefficients which occur in equations (170) and (171). 

A system of only five simultaneous equations for the five 
Fourier coefficients in equation (170) is obtained if first 
equations (163) and (164) are reduced to five equations for 
the five unknowns w, v, u, fi x , and /3«, by means of equations 
(166) to (168). 

For the present, the task is not carried out of obtaining the 
deformation and internal stress Fourier coefficients of 
equations (170) and (171) in terms of the Fourier coefficients 
of the load terms in equation (169). Instead, the axisym- 
metrical case, to which equations (169) to (171) reduce when 
sin md and cos md are interchanged throughout, and then only 
the terms for m = 0 are taken, is treated separately. 

Axisymmetrical deformation of circular cylindrical shell. — 
In equations (163) to (168) set 


d( )/50=o ao/dx = ( )' 

N x e — Qe—Me x = 0 

v=pe=0 me=pe — 0 


( 172 ) 
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and then the following system of equations has to be dealt 
with : 

NJ+p,= 0 1 

Qx — (N el a) +<z=0; 

(173) 

MJ — Q x J r m x ~--Q .. I , . 

<r lm =s—M e /(h+t)a) 

(174) ■ 

0+i>* 

(175) 

N x — vNe =C*u' ) 


Q x —(h-\- 1) Cr c (/3z+ w') 

(176) 

(1 + \)M e - vM x = D*sja E c 1 
M x —vMe = D*\ 3/ ) 

(177) 

The system of equations (173) to (177) may be reduced to 
two simultaneous equations for /3 X and Q x , as follows: First, 
express M x in terms of jS* by means of equation (175) and 
substitute the result in equation (174). From the first of 
equations (177), it follows that 

Me J _|_ x M x + (1 + x)a E ' 

(178) 


and this, introduced into the second of equations (177), gives 


trast with this, the effect of finite G c is to introduce a new term 
into the left sides of these equations. This new term may 
be of appreciable importance, as will be shown. 

Having solved equations (180) and (182), M x and Me are 
obtained from equations (179) and (178), respectively; N> 
follows from equation (173) in the form 

Ne=a(Q x '+q) * (183) 

and w follows from equation (175) in the form 

w=(a/C*) [(1 + X/3 )a QJ + aq+ vf p x dx] (1 84) 

The following examples illustrate the use of equations (178) 
to (184). 

Infinite circular cylindrical shell with periodic load distri- 
bution. — In specialization of equations (169) to (171), set 

q—Qji sin fix s = s M sin /xx ) 

> (185) 

Pz=X>xn cos fix m x =m zit cos p.x) 
w=w„ sin nx u=u„ cos /xx fi x = /3 W cos fix (186) 

Qx= Qm cos nx N X =N X „ sin nx Ne=Neu sin n x ) 

> (187) 

M x =M xlt sin ii x Me=Me„ sin. px ) 

By introducing equations (185) to (187) into equations 
(180) and (182), two simultaneous equations are obtained 
for the amplitudes Q X)i and /3 X(I , as follows: 


M x 


(1+A )D* 
1 — J— X — v 2 




uD*s 

aE c (\+\- v *) 


(179) 


Equation (179) is introduced into the first of equations (174) 
and, restricting attention to shells of uniform section prop- 
erties, there is obtained 


( 1 +X)Z?* n ... i 


vD * fxs^ 


l+X-r 2 

Pxp Jj* [( 1 + g) 


1 + X — v 2 aE 
2 tE f 


a\h T /) G t 


Q. 2 

r Qx»= jj* (m<Z>i+ v PxJE) 


(188) 


(1+X)D* 
1 + X — v 2 


@x Qx — XfL x 


v D*s 

aE c { 1 + X — r 2 ) 


(1 80)* 


To simplify the further discussion, by setting in equation 
(188) m Xfi =p x/1 =0, there is obtained for 0 XI1 and Q x ,, 


To obtain the second of these equations, first, introduce 
into equation (176) the value of w' which follows from 
equation (175), giving 


Qx 

(h + /,) G c 



Ne'-vN/ 


C*(h + t)q' 
12 aE c 



g = L+^ it x 

(1 + X)Z)% 3 

Qxp=— K 


The quantity K is given by 


(189) 


In equation (181), Ne and N x are taken from equation 
(173) and, after slight transformations, there is obtained 


K=\ 1 +-+— — — ^ 

L .3^x 2 (h+t)a 2 G c 


1 — {— X — v 2 
(1 + X)tt 4 


4 1* 


(h- f-Z) 


J' 


(190) 


C*( 1 + 3) e/ (h+t)G c + l3x ~ C7*( 2 ' + ?) (182) * 

Comparing equations (180) and (182) with the correspond- 
ing equations without the effect of transverse shear and nor- 
mal stress deformation, it is seen that the effect of transverse 
normal stress, which is represented by X, merely somewhat 
modifies some of the coefficients of the left sides of the 
corresponding system of equations with E c = 00 . In C on- 


where use has been made of the relation Id equation 

(190) the term X/3 will usually be of little importance. The 
other two variable terms represent the effect of transverse 
shear deformation and of shell curvature, respectively. 
When the radius a is so large that Z 4 /(/i+t) 2 a 2 <d, the shell 
behaves under the action of the given load essentially as a 
plate strip. The effect of transverse shear is important as 
soon as the term (2/ir 2 )(l 2 /a 2 )[t/(h-\-t)](E f /G c ) is not small 
compared with 1. 
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Before evaluating a numerical example the following fur- 
ther formulas which are readily obtained from equations 
(179), (183), and (184) are listed: 


The correction in this case amounts to about 6 percent. 

For a second numerical example, change the above 
moduli ratios to E f /G e = 2000, E f !E c = 1000. This gives 


r 1 

Ne^aq^l —K) 1- 




(191) 


Equations (191) show that in tbis problem not only is the 
deflection increased because of the effect of transverse shear, 
and with that the hoop stress resultant iV# M , but now also an 
effect is found on the bending-moment distribution M XI1 , in 
the opposite sense. The effect of transverse shear is to reduce 
the magnitude of the bending moments in the shell. This 
result is in contrast with what was found for the examples 
which were worked out in the sections on plate analysis and 
circular ring analysis and is therefore of particular sig- 
nificance. 

Equation (191) for w „ may be compared with the corre- 
sponding expression for a simply supported plate strip of 
width l, with sinusoidal load. The result for this case must 
follow from equation (191) in the limit <*> and agree with 
equation (96), which was previously obtained. To compare 
the last of equations (191) with equation (96), the last of 
equations (191) is written in the form 

1 X — v 2 / Z Y 1 , AY 1+A 1 (h+t)tE f 

(1 + X)£>*W w 1 + X — v 2 2 P G c „ 

1 + X-r 2 / X /IV 2 1 E, ^ 
1+ (1 + X)Z>*W a 2 3 + W h+tG c 

~V 1+-K-V 2 /IYC* 

1 (1 + a 2 

(192) 


Equation (192) reduces to the equivalent of equation (96) 
if in it a— > °° . 

From a comparison of equations (192) and (96), it is 
further concluded that the correction due to transverse shear 
is greatest in this case when ct= co, so that, in this case, the 
curvature of the sh ill tends to reduce the additional shear 
deformation, below the value obtained for the simply sup- 
ported plate strip. 

For a first numerical example, take h= 1 inch, <=0.05 
inch, <i — 1 0 inches, 1 = 20 inches, E f /G c = 200, E f /E e = 100, 

i»=l/3, and X=^ 1-05X0.05 iqq — q . 025. The factor K of 
Z luu 

equation (190) becomes 


K— (1+0. 08+38. 6 +54. 5) =0.0106 

instead of K=0.01685. The correction in this case amounts 

0.0180 — 0.0106 „ rT11 
to- o ~ oio6 ~X100~70 percent, thus agam a case is 

found where omission of the effect of transverse shear defor- 
mation would give results which could not be used. How- 
ever, it is noted that the effect of transverse normal stress 
deformation is quite small and may here safely be neglected. 

If the foregoing values of K are introduced into equations 
(191), it is seen that the percentage corrections apply to the 
bending-moment value directly but that for hoop tension 
and radial deflection the corrections are very small indeed. 
In fact, in order that there be appreciable corrections due 
to transverse shear on hoop tension and radial deflection, it 
is necessary that the half wave length of the sinusoidal load 
q be so small that K is at least of magnitude 0.25 or more. 

A case of approximately this kind is obtained if the half 
wave length l is changed from 20 to 10 inches and the moduli 
ratios are again taken as E f IE c — 100, E f IC c — 200. Then, 

K =( 1 +0.008 + 0.965 + 3.41)- 1 = 0.1865 


whereas 


(K)« c =c c = o> = ( 1 + 3.41) _1 = 0.227 

The percentage change of K and therewith of M x is slightly 
more than 19. The percentage change of No and w is about 
4.5. 

The foregoing numerical examples show that the effect 
of transverse shear may be significant in cylindrical sand- 
wich-shell analysis and that moreover its magnitude will 
not in general be predictable by the analysis of an equivalent 
flat-plate or straight-beam problem. 

For the infinite circular cylindrical shell with load q = 
q „ cos nx the essential results are given by equations (190) 
and (191). These results may be extended directly to the 
loading condition 


cos 

Hn = nir/l 


(193) 


By superposition, from equation (191) the following formulas 
are obtained : 


— 8 ', (gJnJ)Kn cos }j. n x 


I^fl + 0.008+4 

- 7T 


400X0.05 4 ( 1 - 0 . 09 ) leo+ooi- 1 

100X1.05 X U + 7T 4 100X1. lj 


= (1 + 0. 008 + 3. 86 + 54. 5)-‘ = 0. 01685 


while without transverse shear and normal stress deformation 


(K)®,-*.— = ( l + 54.5)- 1 = 0.0180 


N»=a’^2q n (l—K n ) cos ix n x 
w=(a 2 /C*)^2q n [l — (l + X/3)If„] cosn n x J 
The values of K n are obtained from the formula 
X, 2 IH E r , 4(1 +X — r 2 ) ' V 


(194) 


#*=[ 1 + 


3 n 2 Tr 2 (h-\-t)a 2 G c (l+\)7r 4 n 4 (h-\-t) 2 a 


■r 


(195) 
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Having the solution for the infinite shell with periodic load 
distribution, it will be only necessary to add to this the gen- 
eral solution of the differential equations without external 
load terms, in order to obtain the complete solution for any 
edge condition of the axisymmetrically stressed circular 
cylindrical shell of finite length. This additional solution 

will now.be obtained. « ' 

Finite circular cylindrical shell acted upon by edge mo- 
ments and forces. — To solve equations (180) and (182) 
with right-hand sides equal to zero, equation (182) is differ- 
entiated twice and (i x " is substituted from equation (180). 
This gives 


<7*0 + 3 )^ 


(h+t)G c 


Qz'~ 


1 + X — v 2 

'(1 + \)D* 


Q x = 0 


or 

Qx [v —2 m 2 Q x " + 4 m 2 4 Qz= 0 (196) 


where 


_1 j C*l 2 
1 cl \ (1 -\-\/3)(h-\-t)G e 


aV 


t 


E, 


(1 — h X/3)(/fc -)— <) G c 


(197a) 


4 / <7*/4 1+X-r 2 1 4 / 1+X-t/ 2 

V E 1 D* (1 +X)(1 -f-X/3) (l+A)(l-(-X/3) 

(197b) 


The auxiliary equation corresponding to equation (196) is 


r 4 — 2m 1 2 r 2 +4m 2 4 = 0 (198a) 

or 

r 2 = m 2 ± Vtoi 1 — 4wi 2 4 (198b) 


The solution of equation (196) occurs in two different forms, 
depending on whether r 2 of equation (198b) is real or not. 
According to equations (197) and (198b), r 2 is complex as 
long as 


or 


m l 4 <(4m2 4 


J_ f tE, J 4[1-, 2 /(l+X)] [ 
a 4 |_(A + 0(l+X/3)flU ^(A + f)V(l+X/3)J 


(199a) 


To clarify this condition, neglect X (which is of very little 
importance here) and equation (199a) then becomes 

(199b) 


When equation (199) holds, a quantity k may be defined by 

k=^j m 2 -\ri V4to 2 2 — m, 4 (200) 

and the four roots of the characteristic equation are k, k, — k, 
and — k, where a bar indicates the taking of conjugates. The 
solution of equation (196) may be written 


Qx=Cie kx -\-Cie^ kx -\-G 2 e kx -\-C2e kx ( 201 ) 


Where equation (199) does not hold, which is the case for 
very small values of GJE f only, all four roots of equation 
(198a) are real and of the form 


k j = -yj mi 2 + V m i 4 — 4 m 2 2 
k 2 = — k\ 

k 3 =^m 2 — V wii 4 — 4 m 2 2 
k i =—k 3 


( 202 ) 


and the solution of equation (196) can be taken in the form 


Q x = A ie klX + A 2 e~ k ' x + A s e k * x + A ie ~ k * x (203) 


Before applying either solution to a specific problem, there 
are noted the following relations which follow from equation 
( 200 ) : 


kk=\k\ 2 —2m 2 2 
k + k=y/2^m 1 2 +2m 2 2 . 


(204) 


Semi-infinite shell acted upon by edge bending moment 
and shear force. — The following boundary conditions hold: 


M x (0) = i + x=^ /(0) = M 0 
Qx(G) — Qo 


(205) 


while for x= » these same quantities vanish. (For the same 
problem without the effect of transverse shear and normal 
stress, see reference 9.) 

Of particular interest in this solution are the values of 
deflection w(0) and change of slope /3 X (0) at the section where 
the loads M 0 and Q 0 are applied. 4 

Taking first the case E f /G c <^2alt for which equation (201) 
applies, it is seen that the conditions at infinity require that 


C 2 = <7 2 =0 (206) 

so that 

Q x =C,e- kx +C,e--* (207) 


The values of (3 X may be obtained by integration from 
equation (180) in the form 


(1+X)Z7* 
1 + x — V 2 


«-“+p 


(208) 


where two constants of integration have been discarded to 
satisfy again the conditions at infinity. 

With equations (207) and (208) there is obtained from the 
boundary conditions (equations (205)) that 


C\+G\=Q 0 

(C i /k) + (C l /ic) = -Mo 


(209) 


< Without transverse shear and normal stress deformation these relations are 

w w Qo ] 

° Mo ] 

in agreement with equations (236) of reference 9, where the homogeneous shell is considered. 
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This determines Ci and Oi in the form 

kQ a +\k\*M 0 

/1_ k-k 

7=; _ hQ a +\k\ 2 M„ 

1 k-k . 


( 210 ) 


Equation (210) is introduced into equation (208) and there 
is obtained as the first of two “influence coefficient” formulas 


l+i j^ fe(°) = ' ~|jp [<2o+(fc + ft)M 0 ] (211) 

The second of these formulas follows from equations (184), 
(207), and (210) in the form 

T ^73^ = -[ | fc| 2 Mo+(fc+*)e°] (212) 

Equations (211) and (212) may be written in more explicit 
form, using equations (204) and (197). The results are 


The effect of transverse shear modifies the deflection due to 
Q 0 and the rotation due to M 0 but not the other two coeffi- 
cients, (2) the effect of transverse normal stress enters all 
four coefficients but only in a minor way, and (3) the reci- 
procity relation that the deflection due to M„ is the same as 
the rotation due to Q 0 is carried over from the theory with- 
out the extra effects. 

For a numerical example the following data are chosen: 
£=0.1 inch, h= 1 inch, a= 10 inches, E f /E c — 100, E f /G c — 200, 

i>= 1/3. This makes X=^ 100 = 0.055, and, from 

equation (197), 


m l = 


/ 01 
“10 V 1-01 


X 200 
018X1. 1 : 


= 0.426 




1-0.09 
100X1.21X1.018 


==0.294 


Then, according to equation (204), 

|fc| 2 =0.173 jfc + F= V2 V0. 182 + 0. 173 = 0. 84 


&( 0 ) 


a 


/(I ~ t~ X/3)(l + X — v 2 ) 

V l+x 



+ 



4 C* 1+X-r 2 
a 2 D* (l+X)(l + X/3) 


(h+t)(i+\/zr °J 


(213) 


and 


w(0) = 


® .« /(I + X/3)(1 +X — 

yC^D* V 1 + x 



o 2 /, ,XW I4C* l + x-r 2 C*la*G c n 

<7* \ ~^3/ V V a 2 D* (l + X)(l+X/3) + (/(,+f)(l+X/3) 

(214) 


Neglecting the generally small effect of finite E c in equa- 
tions (213) and (214), that is, putting X=0 in these equa- 
tions, there may be written instead 


/U0) = 


aVl-*T 0 I /Ti tfa E tl - 

L Qo+ V aW* V + 2 ^ 


w(0) = 


a^l-vK a 2 4/4(1 — / 

“ JC*D* ° C* V a 2 D* V 


t/a E r 


(215)* 
Qo 


Gc 
(216)* 


Equations (215) and (216) contain the noteworthy fact 
that the correction factors for the effect of transverse shear 
are independent of the ratio t/h of face-layer thickness to 
core thickness. The complete formulas of course must and 
do contain the influence of the core thickness h. 

It is further noted that, while equations (211) to (216) 
have been derived for the case that m 1 4 <4m ! ( for which the 
complex solution holds, they are also valid, as is readily 
shown, when 4m 2 4 =£m 1 4 . 

Comparing equations (213) and (214), and (215) and (216) 
with the equations listed in footnote 4 it is seen that: (1) 


while without transverse shear deformation (mi = 0) the 
value of &+&=0.59. According to equations (211) and 
(212), the effect of transverse shear in this case is to in- 
crease the rotation due to the edge moment in the ratio 
0.84/0.59 = 1.42, an effect of 42 percent. The same increase 
is found for the deflection due to the edge shear force. Ro- 
tation due to the shear force and deflection due to the 
moments are practically unchanged. Likewise, the effect 
of transverse normal stress in this case is of negligible 
importance. 

As a further numerical example there is chosen £=0.05 
inch, h= 1 inch, a=20 inches, E f /E c = 1000, Ef/G e = 2000. 

1 1 05X0 05 

This makes X =2 ~ 40 0~ — 1000 = 0.065, and, from equa- 


tion (197), 


m. 




0.05 


1.022X1.05 


2000 = 0.483 


m 2 = 


Vl .05X20 


Vt^" 0 - 2 ' 8 


From equation (204) then 

[A: | 2 — 0.095 k + k=^2 V0. 234 + 0. 095 = 0. 82 


while without transverse shear deformation (?«., = 0) the 
value of &+&=0.44. Thus the effect in this case is to in- 
crease edge rotation due to edge moment and edge deflec- 
tion due to edge shear force in the ratio 0.82/0.44=1.87, an 
effect of 87 percent. 

Infinite circular cylindrical shell acted upon by transverse 
line load. — Calculation is restricted to the determination of 
deflection and bending moment at the section x=0 where the 
line load of intensity 2 Q 0 is assumed to act. The result of 
the foregoing paragraph may be used as follows. Consider 
the infinite shell cut in two parts at the section x=0 and 
assume a bending moment M 0 of such magnitude that the 
slope /L(0) is zero. According to equation (211), this gives 
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M 0 = &== 

k-\-k V 2 TOi 2 +4m 2 2 

and therewith 


(217) 


C*w( 0) __ (k + k) 2 — |fc | 2 „ 2m 1 2 +2m 2 2 

(l+X/3)a- k + k 0 72«ii 2 +4m 2 2 


Qo (218) 


Equations (217) and (218) become, with equations (204) and 
(197), 


M 0 = 


V 2 V 


c* 


1 +x 


Qo 

— 


a 2 D* (l+X)(l+X/3) 


C*/a 2 1 

{h-\-t)G c l+X/3 


(219) 


u;(0) 



V 


C* 1 + X — v 2 
a 2 D* (l+X)(l+X/3) 


(77a 2 1 

(h-\-f)G c l + X/3 



(7* 

a 2 D* 


l+\-v 2 
(l + X)(l + X/3) 


(77a 2 1 

(Ji-\-f)G c 1 + X/3 


( 220 ) 


To give these formulas a somewhat less unwieldy appear- 
ance, the effect of finite E c , that is, X— 0, may again be neg- 
lected, as is permissible in most cases; and there may be 
written 


M — V(^+Qa ~ Qo 

0 2 7 ( 1^77 r 1 t E f 
V i_ 2-v / l-r 2a Qc 


(' 221 )* 


w( 0) 



a 2 V 1 — v 2 a Go Qo 

toJ~(h-\-t)a I 1 ( E r E f 

V + 2 Vl °c 


( 222 )* 


Some numerical examples are as follows. 

Taking ( = 0.1 inch, a = 10 inches, E f !G c — 200, and p=l/3, 
transverse shear deformation reduces M 0 to 1/72.05 times 
the value which holds when G c = ; that is, there is about a 

30-percent reduction in M 0 . At the same time the deflection 
under the line load is 3.05/72.05 = 2.14 times what it is when 
G c = “ ; that is, there is an increase of about 115 percent in 
w( 0) . 

Taking ( = 0.05 inch, a = 20 inches, and E f /G c =200, M 0 is 
decreased by a factor 74/5 = 0.89, while w(0) is increased by 
a factor 1.5/71-25 = 1.34. 

Taking (=0.05 inch, a = 20 inches, and E,IG C — 2000, M 0 
is decreased by a factor 1/73.62 = 0.526, while w (0) is 
increased by a factor 6.25/73.62=3.29. 

Equation (220) for w(0) may be compared with equation 
(116) for the circular plate of radius a. This comparison 
shows that, while for the plate both the ratios f/a and ( h-\-t)/a 
enter into the correction factor, the correction factor for the 
cylindrical shell contains the ratio (/a only; that is, the 
corrections (but not the results) are independent of the ratio 
of face-layer thickness to core thickness in this case of a 
cylindrical shell. * 


SPHERICAL SHELLS 

In conformity with customary usage, the following nota- 
tion is introduced: 


ki=a<t> 

£2 

a t = 1 

£*2 = SHI (f) 

Ri — R 2 = cl 

Nn=N* 

N 22 =N e 

n 12 =n 21 = 

N<t>B Qi = Q<t> 

Q 2 — Qe 

m u =m 4 , 

]\&22 == -Me 

M^\ 2 =Al,p0 

Pl=p* 

II 

e* 

m t =m^ 

m 2 — mo 

Ui—U 

II 


P2=Pe 


X (223) 


Attention is here restricted to problems with rotational 
symmetry and the following relations are used: 


b( )/£>0= 0 

N 4> e=Qe=M 4 ,e=pe=me = v = (3e=0 


(224) 


The differential equations of equilibrium (25) to (28), (31), 
and (34) become, setting 

5( )/&*=d( )/d*=( )' 

(sin <f> Nt)'— cos 4> IVe + sin <j> Q^+a sin <f> p$= 0 (225) 

(sin <(> Qi)' — sin <j> (N^, + Ne) + a sin <f> g = 0 (226) 

(sin <j> M+)' — cos <t> Me — a sin <j> Q^+a sin <j> m^= 0 

(227) 

= 0 (228) 


The stress-strain relations (56), (59) to (61), and (64) to 
(67) become, if there is set in accordance with equation (63) 


Xl X 2 — X 12 — 


2 a 2 “ E e 


(229) 


( 1+ |) W ,-(,-|) N< _0.(^ + f±/A) (230) 

< 231 > 

Q*=Qi+t)G c (a+- (232) 
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(1 + X)M*-(r-X)M»=^(W+-j0 (233) 

(l + X)Ms — ( v — 3<t, cot <£+-gr^ (234) 


There is first given a simple special solution of this system 
of equations and then a generalization is obtained of the two 
simultaneous equations for and (3^ which are fundamental 
in the theoiy of homogeneous isotropic shells. 

Uniform stress distribution in a spherical shell. — Set in 
equations (225) to (234) p^=m^, = 0 and assume that A+ Ne, 
Q 4 , t M#, and Me are independent of 4>. From equation (225) 
it follows that: 

Nt=Ne=N 0 Q* = o (235) 

From equation (226) it follows then that 

iv 0 =! aq (236) 


and from equation (227) it follows that 


II 

£ 

II 

(237) 

Equation (228) gives 


<T{ m = s— 2 M 0 /(h J rt)a 

(238) 


In equations (230) and (231) set u= 0 for reasons of sym- 
metry and obtain 



\ A , <?* 
V ) N *=- 


W 


t(h+t) E f 
6 a E c q 


(239) 


or, with N 0 from equation (236) and X from equation (229), 
(C*/a)w=(l/2)(l — v)aq (240) 


Equation (232) is identically satisfied when /J*=0. Equa- 
tions (233) and (234), in conjunction with equation (237), 
give 


(1 + 2X- 



_s___l t(h-\-t) 2 Ef 
E c 2 a E C S 


or 


M 0 


( h -\- 1) a X 
1 -j- 2X — v 


s 


(241) 


Then, from equation (238), 

(I — *)* 

1 + 2X — v 


(242) 


Equation (242) may be compared with equation (136) for 
the circular ring. 

According to equation (46), there are obtained from equa- 
tions (236) and (241) the following expressions for the stress 
resultants in the outer (“upper”) and inner (“lower”) face 
layers : 


( 1 1 h + t\ 

£ 

II 

© 

f( L\ 

Xs 

V + 2a) 

V4 H 

1+2X — v 


( h+t X 

| N t =a 

(i- 

Xs \ 

V 2a , 

u 

1+2 \ — vJ 

✓ 


Comparison of these results with the corresponding results 
for the circular ring (equations (135)) shows that for given 
values of q and s there is a greater difference between N u and 
Ni in the spherical shell than there is in the circular ring, 
the reason being the relatively larger influence of the s-term 
in equation (243). 

For a specific example, it is again assumed that the radial 
load is applied entirely to the inner face so that g„=0 and, 
according to equations (20) and (22), 



Substitution of equation (244) in equation (243) gives 


0 +&>-( 1 -£ 0 , ¥ 

AT _/i h + L '\ a, J‘ 1+4X — V 

Nl V 2a / 4 1 + 2X — v 


1—v 

1 + 2X — v 

>~ 


(245) 


As a numerical example, taking X= 0.0595, as in the ex- 
ample given in the section entitled “Closed circular ring 
acted upon by uniform radial load,” and r=l/3, it is found 
that the factor in N t which contains the effect of the core 
flexibility is (l + 0.36)/(l + 0.18) = 1.15. Thus, where for the 
circular ring there was a 6-percent stress increase, there now 
is a 15-percent stress increase for the spherical shell. 

Reduction of axisymmetrical problem to two simultaneous 
equations for Q, P and /3*. — The fundamental results of refer- 
ence 10 for homogeneous shells may be readily extended to 
sandwich shells, as follows: 

Equations (225) and (226) are used to express N$ and Ne 
in terms of Qf. 


N 4 , = cot 4> Ei(4>) (246) 

N,= QS+FAt>) (247) 

In equations (246) and (247) the functions F, and F 2 are 
given by 


F,= . % ■ f{q cos 4 >— p# sin <t>) sin <p d<j> (248) 
sin 2 

F 2 = — - - [(sin <(> FiY + apt sin 4 >] (249) 

cos <p 

Next the displacement components u and w are expressed 
in terms of Q$, by means of equations (230), (231), (246), 
and (247). 

Subtraction of equation (231) from equation (230) gives 

— (u r — cot <p u) = (l + v) (N^, — Ne) 
a 

= (1 + v) [— (Q /— cot <f> QY)-\- Fi — F 2 ] 

(250) 


111 1 n iiina 
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Equation (112) is integrated to 


(<7>)u=-(l + r)(<2* + E 3 ) 


where F 3 is given by 


— sin <f> | 


FM — Fifa) 


Equations (251) and (252) are introduced into equation 
(231) and the following expression is obtained for w: 

(C*/a)w= (1 + X/3) (cot <f> Qt+Qt')+F< (253) 

where F 4 is given by 


__ (h+t) JL 
4 12 E, 


+ (l + v) cot $ E 3 +^l +3^) F 2 — 3^ F 


Equations (251) aud (253) are introduced into equation 
(232) for Qj, and the first of the two simultaneous equations 
for <2<#> an d fit is obtained in the form 


(h + 1) G c 


=fif\"fT* [^1 +^(cot 4> Qt~\~ Q/Y + F V + 

■ v ) (Qt + F 3 ) J 


which may be rearranged to read 


Q/' + cot <j> Qy [cof* i+x/3 + ff, (h +0(1+ X/3)_ + 


l+X/3 


fit — FM 


When X = 0 and G c = » and when no external loads are 
present, this equation checks with the first of equations (g) 
on page 469 of reference 9. 

The function F 5 is given by 

v _(1 + v)F i J r F 4 ('oka'i 

l+X/3~ (2 ° 6) 

Introducing the operator 

L = ( )' + cot <t>( )' — cot 2 <£( ) 
equation (255) may finally also be written 


where 


( L — fit — Fs(<t>) 


( 2 t E, X\ 1_ 

1 \h+t G c v ^3j 1 + X 


The second of the two simultaneous equations is obtained 
somewhat more directly as follows: Write equations (233) 
and (234) in the form '• 


D*/a ‘ 

1 — r 2 T 2X(1 T v) 


(l + X)d^. , + (r — X) COt <p fit + 


Me— i_ j/2 ^2X(l-|-r) [ ( ' 1 cofc ^ P'f’ + O' - X)/V + 

(i + r)-|] (260) 

Introduce equations (259)- and (260) into the moment equili- 
brium equation (227) and obtain 

1 - r 2 + 2X(1 + v) [^" + cot <t> fit' ~(™t 2 <I>+y+x) &* + 

(l + r) aQt+amt—0 (261) 

Again, using the operator L, this may be written in the form 

(L- fit- jj* [1 ~ ^+2X(1 + „)] Ffa) (262)* 

The function Fa is given by 

E B =-(l + r)-g---^[l-r 2 +2X(l + r)]m^ (263) 

Equation (262) may be compared with the second of 
equations (g) on page 469 of reference 9. 

Analysis of edge effect for spherical shell. — The special 
case of no distributed surface load and no concentrated load 
at the apex of the shell is obtained by setting 

F t =F t =0 

Following again a known procedure from the theory with- 
out transverse stress deformation, there may be set 

Q,t = ^= fi*=-i== (264) 

ysin 4> ysin 4> 


O' «*' 1 rn t Qi 

Ysin $ z ysin <£ 


Qt"=~^h=-c°t *y=+( f cot 2 *+|)t§4 
ysin <j> V sin 4> k 4 ^/ysin <f> 


+++(!■ 

ysm 4> V 4 


with corresponding formulas for fi<j>' and ficj>". Introduction 
of equation (264) into equations (257) and (262) gives 

C>i / '+(Mi+|-|cot 2 ^Q l -l- I: ^g^ 1 = 0 (266) 

■ cofc2 I 1 “ " 2 + 2X ^ + vMi=0 


Assuming that cot <p 0 is not large compared with unity 5 
and that the effect of the edge loads is restricted to a narrow 
edge zone so that | Qx | | Qx" | , | fix | | fix" | , equations (266) 

and (267) may be simplified to 


«■ ''~^*+TTx73 *=° 


[l-r 2 +2X(l + r)] <2 1 = 0 


(269)* 


4 When cot 6{S> 1 the shell is termed a “shallow” shell which is not considered in what 
follows. 
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Equations (268) and (269) show that the influence of finite 
E c (\^ 0) in the edge-effect problem consists, except in ex- 
treme circumstances, in minor modifications of the results 
for E c — co . The quantity ^ which represents the influence 
of finite G c and which is approximately 


„ _ 2 1 Ef 1 _ 

^ 1 h-\-t G c l+X/3 


( 270 ) 


may, however, in practical cases be large compared with 
unity and not of negligible influence on the results. 

Equations (268) and (269) may be compared with equa- 
tions (180) and (182) for the cylindrical shell. This com- 
parison shows that the influence of finite C c in the edge-effect 
problem is of the same nature for the spherical and cylindri- 
cal shells. Thus, results of the same quantitative nature 
will be obtainable as in the section on cylindrical shells under 
the headings entitled “Finite circular cylindrical shell acted 
upon by edge moments and forces” and “Semi-infinite shell 
acted upon by edge bending moment and shear force.” 

This work is not herein carried further to specific applica- 
tions. It is apparent that such applications may be worked 
out with hardly any more difficulty than when the effect of 
the core deformability is not taken into account. 

CONCLUDING REMARKS 

A system of basic equations has been derived for the analy- 
sis of small-deflection problems of sandwich-type thin shells. 
This system of equations reduces to Love’s theory of thin 
shells when the transverse shear and normal stress deforma- 
bility of the core of the sandwich is of negligible importance. 
The system of basic equations has been applied to a number 
of specific problems from the theory of plates, circular rings, 
circular cylindrical shells, and spherical shells, and it has 
been found that the effects of both transverse shear and trans- 
verse normal stress deformation may be of such magnitude 
that an analysis which disregards them gives values for deflec- 
tions and stresses which are appreciably in error. 

Numerical calculations have been in the nature of sample 
calculations, illustrating both the use of the equations and 
the possible effects of using them. Examples have been 
chosen from the point of view of relative simplicity as well 
as with the thought to illustrate most clearly the conse- 
quences of the extra deformations which have been taken 
into account. It is unavoidable that, in so doing, some of 


the examples may be of little interest for aircraft structural 
analysis and that some problems may not have been analyzed 
which would have well fitted within the contents of this 
report and which at the same time would have been of con- 
siderable practical importance. 

The general analysis has been restricted by the following two 
order-of-magnitude relations: (1)</A<C1 and (2) tEf/hE^l, 
where t is the face-layer thickness, h is the core-layer thick- 
ness, E f is the elastic modulus of the isotropic face-layer 
material, and E c is the elastic modulus in the transverse 
direction of the core-layer material. Therewith it is felt 
that very likely nearly all situations have been covered in 
which the effect of transverse core flexibility is of signifi- 
cant practical importance. It is evident, however, that if 
desired the theory could be extended so as to include cases 
where one or both of these two order-of-magnitude relations 
are not satisfied. The main limitation of the present analy- 
sis is the omission of all finite-deflection and instability 
effects. 

Massachusetts Institute of Technology, 

Cambridge, Mass., May 26, 1947. 
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Absolute coefficients of moment " v 1 

L - M ^ ^ p -.N- y'y - -'y ~ "v - position), 5. (Indicate surfacpby proper subscript.) 


. gbS \P Cm gcS y Gn ~gbS 
(rolling) (pitching) (yawing) 


(-'? Angle of sat j>f .control surface' (relative- to neutral 

• ript.) 


K„ 




4. PROPELLER SYMBOLS 




■ y 


4_Z). ' Diameter-- 

p -- -‘ Geometric pitch ~ ^ . -C- 

p/Dyp’Pitfih. Tatio -S . , . 

y . Inflow velocity 

V, Slipstream velocity . — 

T. Tfirust, absolute ‘Coefficient . 6^=^A=j7 

■> --=r*y;. ■: ■; - - pn 2 D* ; ^ . ,„y, y 

0 - ' Effective helix angle =tany M s -) ' 

»=- ^ 1 . “\V2rrra/ .. 


P .Power, absolute coefficient C P = — 

' y - pn i D s 

■ - ••••' — . . • r- ' ' 1/nVTpt... •;••• .1 

. y ■ C, Speed-power coefficient =vWp^j 

.. . ' ; v ” rj ■ ' Efficiency, - y- y- -U- •' 

ft Revolutions per second, rps s r > < 


d8 < ue , .absolute coefficient Cq= 


\ xx-x. — ;V 


■ -c 




5. NUMEKICAL RELATIONS 


1 hp=76.04 kg-m/s— 550 ft-lb/sec 
1 metric horsepower— 0(P863 hp - 
-;Pmph= 0.4470 mps. , . • *£_>_'- - -r 
p r 1 mps%2.2369 mpfi s - “ . ... 


• 1 lb=0.4530 kg ' 

. I 1 kg=2.2046 lb ' - ‘ 

l mi=l,609,35 m= 5,280 ft 
1 m=3.2808 ft r ; . .: 
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